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Abstract 
Background 


We have previously demonstrated that ivermectin used as prophylaxis for coronavirus disease 2019 (COVID- 
19), irrespective of the regularity, in a strictly controlled citywide program in Southern Brazil (Itajaí, Brazil), 
was associated with reductions in COVID-19 infection, hospitalization, and mortality rates. In this study, 
our objective was to determine if the regular use of ivermectin impacted the level of protection from COVID- 
19 and related outcomes, reinforcing the efficacy of ivermectin through the demonstration of a dose- 
response effect. 


Methods 


This exploratory analysis of a prospective observational study involved a program that used ivermectin at a 
dose of 0.2 mg/kg/day for two consecutive days, every 15 days, for 150 days. Regularity definitions were as 
follows: regular users had 180 mg or more of ivermectin and irregular users had up to 60 mg, in total, 
throughout the program. Comparisons were made between non-users (subjects who did not use ivermectin), 
and regular and irregular users after multivariate adjustments. The full city database was used to calculate 
and compare COVID-19 infection and the risk of dying from COVID-19. The COVID-19 database was used 
and propensity score matching (PSM) was employed for hospitalization and mortality rates. 


Results 


Among 223,128 subjects from the city of Itajaí, 159,560 were 18 years old or up and were not infected by 
COVID-19 until July 7, 2020, from which 45,716 (28.7%) did not use and 113,844 (71.3%) used ivermectin. 
Among ivermectin users, 33,971 (29.8%) used irregularly (up to 60 mg) and 8,325 (7.3%) used regularly 
(more than 180 mg). The remaining 71,548 participants were not included in the analysis. COVID-19 
infection rate was 49% lower for regular users (3.40%) than non-users (6.64%) (risk rate (RR): 0.51; 95% CI: 
0.45-0.58; p < 0.0001), and 25% lower than irregular users (4.54%) (RR: 0.75; 95% CI: 0.66-0.85; p < 0.0001). 
The infection rate was 32% lower for irregular users than non-users (RR: 0.68; 95% CI: 0.64-0.73; p < 0.0001). 
Among COVID-19 participants, regular users were older and had a higher prevalence of type 2 diabetes and 
hypertension than irregular and non-users. After PSM, the matched analysis contained 283 subjects in each 
group of non-users and regular users, between regular users and irregular users, and 1,542 subjects between 
non-users and irregular users. The hospitalization rate was reduced by 100% in regular users compared to 
both irregular users and non-users (p < 0.0001), and by 29% among irregular users compared to non-users 
(RR: 0.781; 95% CI: 0.49-1.05; p = 0.099). Mortality rate was 92% lower in regular users than non-users (RR: 
0.08; 95% CI: 0.02-0.35; p = 0.0008) and 84% lower than irregular users (RR: 0.16; 95% CI: 0.04-0.71; p = 
0.016), while irregular users had a 37% lower mortality rate reduction than non-users (RR: 0.67; 95% 

CI: 0.40-0.99; p = 0.049). Risk of dying from COVID-19 was 86% lower among regular users than non-users 
(RR: 0.14; 95% CI: 0.03-0.57; p = 0.006), and 72% lower than irregular users (RR: 0.28; 95% CI: 0.07-1.18; p = 
0.083), while irregular users had a 51% reduction compared to non-users (RR: 0.49; 95% CI: 0.32-0.76; p = 
0.001). 


Conclusion 


Non-use of ivermectin was associated with a 12.5-fold increase in mortality rate and a seven-fold increased 
risk of dying from COVID-19 compared to the regular use of ivermectin. This dose-response efficacy 
reinforces the prophylactic effects of ivermectin against COVID-19. 
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Introduction 


Ivermectin has been proposed as potential prophylaxis and therapy for coronavirus disease 2019 (COVID- 
19) due to its previously reported anti-viral [1-4], metabolic [5-10], and anti-inflammatory [11-19] actions, 
with strong plausibility [20,21] and positive in vitro, in vivo, and epidemiological findings [22-24] in 
preliminary studies. 


Between July and December 2020, a citywide program in Itajai, in the state of Santa Catarina, Southern 
Brazil, offered a voluntary, medically prescribed program of ivermectin as prophylaxis for COVID-19. This 
was based on the extensive, well-established safety profile and known absence of risks with long-term use of 
ivermectin, and the lack of therapeutic and preventive alternative options in 2020. 


The systematically collected data within this program demonstrated that ivermectin used as prophylaxis for 
COVID-19 improved COVID-19 related-outcomes. The use of ivermectin led to a 44% reduction in infection 
rate, a 56% reduction in hospitalization rate, and a 68% reduction in mortality rate by using propensity score 
matching (PSM) to balance the study groups [25]. 


These conclusions were based on an analog evaluation of the intent-to-treat (ITT) analysis of randomized 
clinical trials (RCTs). All participants of the program were included for analysis, irrespective of regularity or 
the total amount of ivermectin taken. Among participants of the ivermectin use (regular and irregular) as 
prophylaxis for the COVID-19 program, it was unknown if regular ivermectin use would lead to a more 
substantial reduction in COVID-19 infection rate and related outcomes than irregular use. 


In this study, an evaluation was done with participants that used ivermectin prophylactically for COVID-19, 
to determine if regular use compared to irregular use impacted the degree of reduction in COVID-19 
infection, hospitalization, and mortality rates. Regular and irregular ivermectin users were also compared to 
non-users to evaluate evidence of a dose-response pattern of efficacy. 


Materials And Methods 
Study population 


A thorough description of the program, study population, and protocol was described elsewhere [25]. This 
was a medically based, observational, and prospective study that involved the voluntary use of ivermectin as 
prophylaxis for COVID-19 in the city of Itajai, Santa Catarina, Brazil. It was a citywide program conducted 
between July 7 and December 2, 2020. Data were collected prospectively and systematically, as was the 
mandatory reporting of all events. 


The study design, institutional review board (IRB) approval, and data analysis were done upon completion of 
the program. The study of the COVID-19 cases reported in the city of Itajai (n = 9,956, including cases that 
occurred before July 7, 2020, as a comparison) was approved by the National Research Ethics Council 
(CONEP) (approval number: 4.821.082; protocol (CAAE) number: 47124221.2.0000.5485). 


Study procedures and data collection 


Voluntary prophylactic use of ivermectin was offered as an option to patients during medical visits in a 
provisional outpatient clinic at the Convention Center and in secondary outpatient clinics at local health 
centers in the city of Itajai, as part of the Universal Health System (SUS). During medical visits, patient data, 
including medical history, comorbidities, previous diseases, medications, and physical signs (body weight, 
height, body mass index, systolic and diastolic blood pressure, and heart rate), were recorded in the SUS- 
based system. Ivermectin was then optionally prescribed in a dose of 0.2 mg/kg/day for two consecutive days, 
every 15 days to participants who presented without symptoms of COVID-19 or any contradictions to 
ivermectin. 


During the study, subjects who became infected with COVID-19 and diagnosed with a positive reverse 
transcription-polymerase chain reaction (RT-PCR) for SARS-CoV-2 were documented and medically followed 
up. Data on hospitalizations and deaths due to COVID-19 were also systematically registered. 


In this analysis, all residents from the city of Itajai were considered. This included participants in the 
program that used and did not use ivermectin prophylactically. Registry data were analyzed for all 
participants included in the sample. Subjects with a positive diagnosis of COVID-19 before July 7, 2020, 
when the program was initiated, and those below 18 years old were excluded from the analysis. 


The 223,128 residents from Itajaí included 114,568 participants aged 18 years and above, who used 
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ivermectin prophylactically and 45,716 who did not use ivermectin, throughout the citywide program. 
Among these participants, 113,844 were not infected prior to July 7, 2020. This program also included 8,352 
subjects aged 18 years and above from other cities that participated in the program, although not included 
in the present analysis. 


While ivermectin non-users remained unchanged from the first analysis [25], ivermectin users were divided 
according to the accumulated dose of ivermectin taken. 


The analysis focused on data for participants that used up to 60 mg (10 tablets) of ivermectin and those that 
used more than 180 mg (more than 30 tablets). Grouping the users in this manner represented a higher 
certainty of regularity and irregularity, respectively. These groups were compared to non-users in a three- 
group comparison analysis. 


The three two-group matching of ivermectin, i.e., (1) non-users and regular users, (2) non-users and 
irregular users, and (3) regular users and irregular users, were balanced and matched using PSM with the 
following variables: age, sex, history of smoking, myocardial infarction (MI), stroke, hypertension, type 2 
diabetes (T2D), cardiovascular diseases (CVD), cancer (any type), asthma, chronic obstructive pulmonary 
disease (COPD), and other pulmonary diseases. 


Because the accuracy of the reports was guaranteed for Itajaí residents only, all calculations and rates were 
based on the participants from the city. The database used for the calculation of COVID-19 infection rate 
and for risk of dying from COVID-19 was the entire city of Itajai and then calculated among ivermectin 
regular users, irregular users, and non-users of participants from Itajai. Analyses were performed before and 
after adjustment for multiple variables. 


Hospitalization and mortality rates were analyzed for all participants reported with a positive COVID-19 
diagnosis from Itajaí. Reports of all COVID-19 deaths were mandatory, while hospitalization rates were 
based on the data from the local public hospital only, which may justify potential discrepancies between 
hospitalization and mortality rates. We calculated hospitalization and hospitalization rates before matching 
and after PSM of groups, followed by a multivariate-adjusted analysis of the residual differences (double- 
adjusted model). 


In Supplement Appendix 1, pre-matched comparisons of hospitalization and mortality rates are 
provided. Figure / illustrates the locations of each analysis performed in this study. Datasets are publicly 
available at https://osf.io/uxhaf/. 


2022 Kerr et al. Cureus 14(8): e28624. DOI 10.7759/cureus.28624 3 of 37 


Cureus 


CATEGORIES ACCORDING TO LEVEL OF IVERMECTIN USE 


NON -USERS REGULAR USERS IRREGULAR USERS 
(180mg or more) (up to 60mg) 
45,716 8,325 33,971 


3-GROUP PRE-MATCHING — TABLE 1 


Baseline 
characteristics 


3 x 2-GROUP POST-MATCHED — TABLE 2 


DATABASE 
ALL COVID-19 CASES 


3-GROUP PRE-MATCHING — FIGURE 2 


COVID-19 
infection 
rate 


DATABASE 
CITY OF ITAJAÍ 


3-GROUP PRE-MATCHING* — 
TABLE 9/1S AND 10/2S | FIGURE 7/1S 


rate 


COVID-19 


3 x 2-GROUP POST-MATCHED — 
TABLES 3, 4 AND 5 | FIGURE 3 


DATABASE 
ALL COVID-19 CASES 


3-GROUP PRE-MATCHING* — 
TABLE 11/3S AND 12/45 | FIGURE 8/2S 


rate 


3 x 2-GROUP POST-MATCHED — 
TABLE 6, 7 AND 8 | FIGURE 4 


COVID-19 
mortality hospitalization 


DATABASE 
ALL COVID-19 CASES 


3-GROUP PRE-MATCHING — FIGURE 5 


Risk of 


dying 
from COVID-19 


DATABASE 
CITY OF ITAJAÍ 


*Supplement Appendix 1 


FIGURE 1: Illustrative study guide 


Statistical analysis 


The risk of hospitalizations and deaths were calculated for all three groups before matching and for each of 
the three two-group combinations that were propensity score matched. 


Comparisons between groups for hospitalization and mortality rates were calculated using chi-square before 
adjusting for variables and after multivariate adjustments. The generalized linear mixed model was 
employed, assuming the binomial distribution for the residues and included the fixed classificatory effects 
for each of the variables. While there were no missing data, as per the system, illogical data, restricted to age, 
were corrected individually, although some may remain due to the exceptional amount of data gathered. Age 
between 100 and 115 years old was considered presumably illogical, rechecked, and corrected when needed. 
Age above 115 years old was considered obligatorily illogical and corrected accordingly. Statistical analysis 
software SAS/STAT (SAS Institute Inc., Cary, North Carolina, USA) was used for the present study. 
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Results 


There were 159,560 participants aged 18 years and above not infected with COVID-19 prior to July 7, 2020, 
from the city of Itajai, Brazil. Among them, 45,716 (28.7%) did not use ivermectin and 113,844 (71.3%) used 
ivermectin prophylactically. Of the 113,844 participants, 8,325 (7.3%) subjects used ivermectin regularly and 
33,971 (29.8%) used ivermectin irregularly. In total, 88,012 subjects were included in the present analysis. 
The 71,548 (62.8%) remaining participants used intermediate doses between 60 mg and 180 mg and were not 
included in this analysis. 


Before matching, a total of 7,228 subjects from the city of Itajai were infected with COVID-19 between July 7 
and December 2, 2020. Of these, 3,034 (42.0%) did not use ivermectin prophylactically, 283 (3.9%) used 
ivermectin regularly, 1,542 (21.3%) used ivermectin irregularly, and 2,369 (32.8%) used intermediate doses 
of ivermectin. Comparisons between ivermectin non-users, regular users, and irregular users are described 


in Table /. 

Characteristics yea? (n= sy users (n = ee yal users (n= Hees lial the three 
Mean (SD) 39.8 + 14.2 47.0 + 14.2 41.0 + 14.5 
Age <0.0001 
<30 years 844 (27.8%) 39 (13.8%) 397 (25.7%) 

30-50 years 1,582 (52.2%) 131 (46.3%) 775 (50.3%) 

>50 years 608 (20.0%) 113 (39.9%) 370 (24.0%) 

Sex 0.19 
Female 1,624 (53.5%) 141 (49.8%) 853 (55.3%) 

Male 1,410 (46.5%) 142 (50.2%) 689 (44.7%) 

Race 0:055 
Afro-Brazilian 100 (3.3%) 4 (1.4%) 37 (2.4%) 

Mixed 682 (22.5%) 58 (20.5%) 373 (24.2%) 

Caucasian 2,192 (72.5%) 221 (78.1%) 1,102 (71.5%) 
Asian-Brazilian 60 (51.7%) 0 (0.0%) 30 (1.9%) 
Type 2 diabetes 0.33 
Yes 63 (2.1%) 9 (3.2%) 40 (2.6%) 

No 2,971 (97.9%) 274 (96.8%) 1,502 (97.4%) 

Hypertension 0:15 
Yes 166 (5.5%) 23 (8.1%) 96 (6.2%) 

No 2,868 (94.5%) 260 (91.9%) 1,446 (93.8%) 
Asthma 0.47 
Yes 6 (0.2%) 0 (0.0%) 6 (0.4%) 

No 3,028 (99.8%) 283 (100.0%) 1,536 (99.6%) 

COPD 0.42 
Yes 6 (0.2%) 1 (0.4%) 1 (0.1%) 

No 3,028 (99.8%) 282 (99.6%) 1,541 (99.9%) 


Other respiratory 


, 0.78 
diseases 
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eee 5 ER (n= ROTEN users (n = OT users (n = P-value (between the three 
No 088 009 8%) FE o9 6%) 13609 8%) groups) 
ee om 
Yes 15 (0.5%) 2 (0.7%) 16 (1.0%) 

No 3,019 (99.5%) 281 (99.3%) 1,526 (99.0%) 

Cancer 0.73 

Yes 12 (0.4%) 2 (0.7%) 6 (0.4%) 

No 3,022 (99.6%) 281 (99.3%) 1,536 (99.6%) 

History of smoking 0.81 

Yes 47 (1.5%) 3 (1.1%) 23 (1.5%) 

No 2,987 (98.5%) 280 (98.9%) 1,519 (98.5%) 

History of stroke 0.71 

Yes 10 (0.3%) 1 (0.4%) 3 (0.2%) 

No 3,024 (99.7%) 282 (99.6%) 1,539 (99.8%) 

History of MI 0.64 

Yes 4 (0.1%) 0 (0.0%) 3 (0.2%) 

No 3,030 (99,9%) 283 (100.0%) 1,539 (99.8%) 


TABLE 1: Pre-matched baseline characteristics of ivermectin non-users, regular users, and 
irregular users 


SD = standard deviation; COPD = chronic obstructive pulmonary disease; MI = myocardial infarction. 


Baseline characteristics 


Table 1 describes the baseline characteristics of the groups of ivermectin non-users (n = 3,034), regular users 
(n = 283), and irregular users (n = 1,542) before matching groups. Age was significantly different across 
groups for levels of ivermectin use (p < 0.0001). Ivermectin regular users had a higher percentage of subjects 
above 50 years old (39.9%) than irregular users (24.0%) and non-users (20.0%). There were fewer subjects 
below 30 years old among regular users (13.8%) than among irregular users (25.7%) and non-users (27.8%). 
All other baseline characteristics were numerical but not statistically different. There were slightly more 
males among regular users (50.2%) than irregular users (44.7%) and non-users (46.5%) (p = 0.19). The 
percentage of participants with T2D was numerically higher among regular users (3.2%) than irregular users 
(2.6%) and non-users (2.1%) (p = 0.33). Hypertension was more prevalent in regular users (8.1%) than 
irregular users (6.2%) and non-users (5.5%) (p = 0.15). 


Table 2 describes the baseline characteristics of ivermectin non-users paired with regular users and non- 
users paired with irregular users. After balancing and matching each of the three combinations of two 
groups (non-users and regular users, non-users and irregular users, and regular and irregular users), there 
were 283 subjects in each group (n = 566) between non-users and regular users and between irregular and 
regular users, and 1,542 (n = 3,084) between non-user and irregular users, with similar baseline 


characteristics. 
Non-users paired with regular Non-users paired with irregular Regular users paired with irregular 
ivermectin users ivermectin users ivermectin users 


Non-users (n Regular users (n Non-users (n Irregular users (n Regular users (n Irregular users (n 


Variabl 

rons = 283) = 283) = 1,542) = 1,542) = 283) = 283) 
Age 
Mean (SD) 4164148  47.0+142 40.34144 4104145 47.0 + 14.2 43.8 + 16.0 
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Age 

<30 years 
30-50 years 
>50 years 

Sex 

Female 

Male 

Race 
Afro-Brazilian 
Mixed 
Caucasian 
Asian-Brazilian 
Type 2 diabetes 
Yes 

No 
Hypertension 
Yes 

No 


Asthma 


Other respiratory 
diseases 


Yes 
No 


Cardiovascular 
diseases 


Yes 
No 


Cancer 


No 


History of 
smoking 


Yes 
No 
History of stroke 


Yes 


Non-users paired with regular 


ivermectin us 
63 (22.3%) 


152 (53.7%) 


68 (24.0%) 


156 (55.1%) 


127 (44.9%) 


9 (3.2%) 
58 (20.5%) 
213 (75.3%) 


3 (1.1%) 


10 (3.5%) 


273 (96.5%) 


21 (7.4%) 


262 (92.6%) 


0 


283 (100.0%) 


0 (0.0%) 


283 (100.0%) 


1 (0.3%) 


282 (99.7%) 


1 (0.3%) 


282 (99.7%) 


2 (0.7%) 


281 (99.3%) 


2 (0.7%) 


281 (99.3%) 


1 (0.3%) 


ers 
39 (13.8%) 
131 (46.3%) 


113 (39.9%) 


141 (49.8%) 


142 (50.2%) 


4 (1.4%) 
58 (20.5%) 
221 (78.1%) 


0 (0.0%) 


9 (3.2%) 


274 (96.8%) 


23 (8.1%) 


260 (91.9%) 


0 


283 (100.0%) 


1 (0.3%) 


282 (99.7%) 


1 (0.3%) 


282 (99.7%) 


2 (0.7%) 


281 (99.3%) 


2 (0.7%) 


281 (99.3%) 


3 (1.1%) 


286 (98.1%) 


1 (0.3%) 
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z 
1D 
© 
3 
D 
Q 
A 
5 


i us 
410 (26.6%) 


808 (52.4%) 


324 (21.0%) 


846 (54.9%) 


696 (45.1%) 


45 (2.9%) 
351 (22.8%) 
1,114 (72.2%) 


32 (2.1%) 


37 (2.4%) 


1,505 (97.6%) 


86 (5.6%) 


1,456 (94.4%) 


6 (0.4%) 


1,536 (99.6%) 


1 (0.1%) 


1,541 (99.9%) 


3 (0.2%) 


1,539 (99.8%) 


9 (0.6%) 


1,533 (99.4%) 


6 (0.4%) 


1,536 (99.6%) 


21 (1.4%) 


1,521 (98.6%) 


2 (0.1%) 


ers 
397 (25.7%) 
775 (50.3%) 


370 (24.0%) 


853 (55.3%) 


689 (44.7%) 


37 (2.4%) 
373 (24.2%) 
1,102 (71.5%) 


30 (2.0%) 


40 (2.6%) 


1,502 (97.4%) 


96 (6.2%) 


1,446 (93.8%) 


6 (0.4%) 


1,536 (99.6%) 


1 (0.1%) 


1,541 (99.9%) 


3 (0.2%) 


1,539 (99.8%) 


16 (1.0%) 


1,526 (99.0%) 


6 (0.4%) 


1,536 (99.6%) 


23 (1.5%) 


1,519 (98.5%) 


3 (0.2%) 


Regular users paired with irregular 
vermectin users 


39 (13.8%) 
131 (46.3%) 


113 (39.9%) 


141 (49.8%) 


142 (50.2%) 


4 (1.4%) 
58 (20.5%) 
221 (78.1%) 


0 


9 (3.2%) 


274 (96.8%) 


23 (8.1%) 


260 (91.9%) 


0 


283 (100.0%) 


1 (0.3%) 


282 (99.7%) 


1 (0.3%) 


282 (99.7%) 


2 (0.7%) 


281 (99.3%) 


2 (0.7%) 


281 (99.3%) 


3 (1.1%) 


280 (98.9%) 


1 (0.3%) 


60 (21.2%) 
132 (46.4%) 


91 (32.2%) 


155 (54.8%) 


128 (45.2%) 


5 (1.8%) 
68 (24.0%) 
209 (73.9%) 


1 (0.3%) 


10 (3.5%) 


273 (96.5%) 


20 (7.1%) 


263 (92.9%) 


0 


283 (100.0%) 


0 


283 (100.0%) 


0 


283 (100.0%) 


5 (1.8%) 


278 (98.2%) 


2 (0.7%) 


281 (99.3%) 


1 (0.3%) 


282 (99.7%) 


1 (0.3%) 
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No VAENE pairéd with regular INSP AREA aired WH fregtlar *Reditar Geers paired with irregular 
ivermectin users ivermectin users ivermectin users 

History of MI 

Yes 0 0 1 (0.1%) 3 (0.2%) 0 0 

No 283 (100.0%) 283 (100.0%) 1,541 (99.9%) 1,539 (99.8%) 283 (100.0%) 283 (100.0%) 


TABLE 2: Baseline characteristics of the prophylactic study after propensity score matching 
between non-users and regular users, non-users and irregular users, and irregular users and 
regular users 


SD = standard deviation; COPD = chronic obstructive pulmonary disease; MI = myocardial infarction. 


Impact of ivermectin on infection rates in non-users, regular users, and 
irregular users 


Figure 2 illustrates infection rates for ivermectin non-users, regular users, and irregular users, during the 
overall, first, and second half of the program. In the program, the infection rate among ivermectin non-users 
was 6.64% (3,034/45,716 infections). Ivermectin regular users had a reduction of 49% in infection rate 
compared to non-users (283/8,325 cases; 3.40% infection rate; risk ratio (RR): 0.51; 95% CI: 0.45 - 0.58; p< 
0.0001). Irregular ivermectin users had a 32% lower infection rate than non-users (1,542/33,971; 4.54% 
infection rate; RR: 0.68; 95% CI: 0.64 - 0.73); p < 0.0001). Ivermectin regular users had a 25% lower infection 
rate than irregular users (RR versus sporadic users: 0.75; 95% CI: 0.66 - 0.85; p < 0.0001). 
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Citywide program of COVID-19 prophylaxis with ivermectin (Itajaí, Brazil) 
Regular and irregular users 


NON | REGULAR IRREGULAR 
IVERMECTIN USE| {IVERMECTIN USE| | IVERMECTIN USE 
45,716 | 8,325 33,971 
OVERALL 
3,034 283 1,542 
__ cases cases cases 
6.64% 3.40% 4.54% REEE 
Pe, AE panei p-value 
infection rate infection rate infection rate 
Group 
got 0.68 | compared _ 
(0.45 — 0.58) (0.64 - 0.73) à m 
< 0.0001 < 0.0001 RR = risk ratio; 
vs non-use Enone CI = confidence interval 
0.75 
(0.66 — 0.85) 
< 0.0001 
vs irregular use 
FIRST HALF SECOND HALF 
UUL7-SEP 19) (SEP 20- DEC 2) 
1,422 | 121 908 1,612 162 634 
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< 0.0001 ; 0.63 vs 
vs irregular use irregular use 


FIGURE 2: Impact of ivermectin use on infection rates during the 
overall, first half, and second half of the program in non-users, regular 
users, and irregular users 


In the first half of the program, between July 7 and September 19, 2020, infection rate was 3.11% (1,422 
cases) among ivermectin non-users and 1.45% (121 cases) among ivermectin regular users; a 53% reduction 
compared to non-users (RR: 0.47; 95% CI: 0.39 - 0.56; p < 0.0001). Infection rate was 2.67% (908 cases) 
among ivermectin irregular users, showing a 14% reduction compared to non-users (RR: 0.86; 95% CI: 0.79 - 
0.93; p = 0.0003). Regular users had 46% lower infection rate than irregular users (RR: 0.54; 95% CI: 0.45 - 
0.66; p < 0.0001). 


In the second half of the program, between September 20 and December 2, 2020, infection rate was 3.53% 
(1,612 cases) among ivermectin non-users and 1.95% (162 cases) among ivermectin regular users; a 45% 
reduction compared to non-users (RR: 0.55; 95% CI: 0.47 - 0.65; p < 0.0001). Infection rate was 1.87% among 
ivermectin irregular users (634 cases), showing a 47% reduction in infection rate compared to non-users (RR: 
0.53; 95% CI: 0.48 - 0.58; p < 0.0001). Regular users had a similar infection rate to irregular users during the 
second half of the program (RR: 1.04; 95% CI: 0.88 - 1.24; p = 0.63). 


Hospitalization rates for ivermectin non-users, regular users, and 
irregular users 


Supplement Appendix 1A shows hospitalization rates before matching. Tables 3-5 show hospitalization rates 
and unadjusted and multivariate-adjusted values for each of the three two-group comparisons after 
balancing and matching. 
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Cureus 


Propensity score matched Ivermectin . E Unadjusted hospital Multivariate adjusted 
non-users and regular non-users (n Regular ivermectin risk ratio (95% Cl) and hospital risk ratio (95% Cl) 
users = 283) users (n = 283) p-value (p) and p-value (p) 

Overall 13/283 (4.6%) 0/283 (0.0%) 0.04 (0.002 — 0.60) (0.02) 0.00 (n/a) (<0.0001) 

Age 

<30 y/o 0/63 (0.0%) 0/39 (0.0%) 1.61 (0.03 — 82.7) (0.81) 1.00 (n/a) (1.00) 

30-50 y/o 3/152 (2.0%) 0/131 (0.0%) 0.16 (0.01 — 3.17) (0.23) n/a 

>50 y/o 10/68 (14.7%) 0/113 (0.0%) 0.02 (0.001 — 0.43) (0.011) n/a (n/a) (<0.001) 

Sex 

Female 7/156 (4.5%) 0/141 (0.0%) 0.07 (0.004 — 1.24) (0.07) n/a (n/a) (<0.001) 

Male 6/127 (4.7%) 0/142 (0.0%) 0.07 (0.004 — 1.18) (0.064) n/a (n/a) (<0.001) 

Race 

Afro-Brazilian 0/9 (0.0%) 0/4 (0.0%) 2/11 (0.04 — 124.5) (0.72) 1.00 (n/a) (1.00) 

Mixed 3/58 (5.2%) 0/58 (0.0%) 0.14 (0.01 — 2.7) (0.19) n/a 

Caucasian 10/213 (4.7%) 0/221 (0.0%) 0.14 (0.01 — 2.68) (0.19) n/a (n/a) (<0.001) 

Asian-Brazilian 0/3 (0.0%) 0/0 n/a n/a 

Type 2 diabetes 

Yes 3/10 (30.0%) 0/9 (0.0%) 0.11 (0.005 — 2.54) (0.17) n/a 

No 10/273 (3.7%) 0/274 (0.0%) 0.05 (0.003 — 0.78) (0.033) n/a 

Hypertension 

Yes 5/21 (23.8%) 0/23 (0.0%) 0.06 (0.003 — 1.24) (0.069) n/a 

No 8/262 (3.1%) 0/260 (0.0%) 0.06 (0.003 — 1.00) (0.05) n/a 

Asthma 

Yes 0/0 0/0 n/a n/a 

No 13/283 (4.6%) 0/283 (0.0%) 0.04 (0.002 — 0.60) (0.02) 0.00 (0.00 — 0.00) (<0.001) 

COPD n/a 

Yes 0/0 0/1 (0.0%) 0.50 (0.04 — 7.10) (0.61) n/a 

No 13/283 (4.6%) 0/282 (0.0%) 0.04 (0.002 — 0.60) (0.021) 0.00 (0.00 — 0.00) (<0.001) 


Other respiratory diseases 
Yes 0/1 (0.0%) 0/1 (0.0%) 1.00 (0.01 — 92.4) (1.00) n/a 
No 13/282 (4.6%) 0/282 (0.0%) 0.04 (0.002 — 0.60) (0.021) 0.00 (0.00 — 0.00) (<0.001) 


Cardiovascular diseases 


Yes 0/1 (0.0%) 0/2 (0.0%) 0.33 (0.02 — 5.33) (0.44) 1.00 (n/a) (1.00) 

No 13/282 (4.6%) 0/281 (0.0%) 0.04 (0.002 — 0.60) (0.021) 0.00 (0.00 — 0.00) (<0.001) 
Cancer 

Yes 1/2 (50.0%) 0/2 (0.0%) 0.20 (0.01 — 8.83) (0.40) n/a 

No 12/281 (4.3%) 0/281 (0.0%) 0.04 (0.002 — 0.65) (0.024) 0.00 (0.00 — 0.00) (<0.001) 
History of smoking 

Yes 0/2 (0.0%) 0/3 (0.0%) 0.71 (0.01 — 49.7) (0.88) 1.00 (n/a) (1.00) 

No 13/281 (4.6%) 0/280 (0.0%) 0.04 (0.002 — 0.60) (0.021) 0.00 (0.00 — 0.00) (<0.001) 
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Hstopensittycicore matched 
non-users and regular 

yers 

No 

History of MI 


Yes 


No 


TABLE 3: Hospitalization rates in the three two-group comparisons after balancing and matching 


Ivermectin 
non-users (n 


OT128Ph%) 


13/282 (4.6%) 


0/0 


13/283 (4.6%) 


Regular ivermectin 


HEH) 283) 


0/282 (0.0%) 


0/0 


0/283 (0.0%) 


Unadjusted hospital 
risk ratio (95% Cl) and 


10434884P) 92.4) (1.00) 


0.04 (0.002 — 0.60) (0.021) 


n/a 


0.04 (0.002 — 0.60) (0.02) 


the groups of non-users and regular ivermectin users 


Multivariate adjusted 
hospital risk ratio (95% Cl) 
rand p-value (p) 


0.00 (0.00 — 0.00) (<0.001) 


n/a 


0.00 (0.00 — 0.00) (<0.001) 


MI = myocardial infarction; COPD = chronic obstructive pulmonary disease; y/o = years old; CI = confidence interval; n/a = not applicable. 


Propensity score matched 
non-users and irregular 
users 


Overall 

Age 

<30 y/o 

30-50 y/o 

>50 y/o 

Sex 

Female 

Male 

Race 
Afro-Brazilian 
Mixed 
Caucasian 
Asian-Brazilian 
Type 2 diabetes 
Yes 

No 
Hypertension 
Yes 

No 


Asthma 


Other respiratory diseases 


Ivermectin 
non-users (n 
= 1,542) 


47/1,542 (3.0%) 


0/410 (0.0%) 
4/808 (0.5%) 


43/324 (13.3%) 


24/846 (2.8%) 


23/696 (3.3%) 


2/45 (4.4%) 
9/351 (2.6%) 
36/1,114 (3.2%) 


0/32 (0.0%) 


6/37 (16.2%) 


41/1,505 (2.7%) 


13/86 (15.1%) 


34/1,456 (2.3%) 


0/6 (0.0%) 


47/1,536 (3.1%) 


0/1 (0.0%) 


47/1,541 (3.0%) 
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Irregular 
ivermectin 


users (n = 1,542) 


38/1,542 (2.5%) 


1/397 (0.3%) 
71775 (0.9%) 


30/370 (8.1%) 


17/853 (2.0%) 


21/689 (3.0%) 


0/37 (0.0%) 
11/373 (3.0%) 
26/1,102 (2.4%) 


1/30 (3.3%) 


3/40 (7.5%) 


35/1,502 (2.3%) 


9/96 (9.4%) 


29/1,446 (2.0%) 


1/6 (16.7%) 


37/1,536 (2.4%) 


0/1 (0.0%) 


38/1,541 (2.5%) 


Unadjusted hospital risk 
ratio (95% Cl) and p- 
value (p) 


0.80 (0.52 — 1.24) (0.32) 


3.11 (0.13 — 76.5) (0.49) 
1.83 (0.53 — 6.28) (0.34) 


0.58 (0.35 — 0.94) (0.028) 


0.70 (0.37 — 1.31) (0.26) 


0.92 (0.50 — 1.68) (0.79) 


0.23 (0.01 — 4.99) (0.35) 
1.15 (0.47 — 2.82) (0.75) 
0.73 (0.44 — 1.21) (0.22) 


3.31 (0.13 — 84.3) (0.47) 


0.42 (0.097 — 1.81) (0.24) 


0.85 (0.54 — 1.35) (0.49) 


0.58 (0.24 — 1.44) (0.24) 


0.86 (0.52 — 1.41) (0.54) 


3.55 (0.12 — 105.8) (0.47) 


0.78 (0.51 — 1.21) (0.27) 


1.00 (0.01 — 92.4) (1.00) 


0.80 (0.52 — 1.24) (0.32) 


Multivariate adjusted hospital 
risk ratio (95% Cl) and p-value 


(p) 


0.71 (0.49 — 1.05) (0.099) 


n/a (0.98) 
1.82 (0.54 — 6.21) (0.34) 


0.61 (0.39 — 0.95) (0.029) 


0.69 (0.38 — 1.26) (0.23) 


0.71 (0.40 — 1.24) (0.22) 


n/a (0.98) 
1.02 (0.44 — 2.34) (0.97) 
0.65 (0.40 — 1.05) (0.078) 


n/a (0.98) 


0.51 (0.14 — 1.85) (0.31) 


0.75 (0.49 — 1.15) (0.19) 


0.59 (0.27 — 1.31) (0.20) 


0.75 (0.47 — 1.23) (0.26) 


n/a (0.58) 


0.70 (0.46 — 1.05) (0.087) 


1.00 (n/a) (1.00) 


0.71 (0.48 — 1.07) (0.11) 
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Yes ; 
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non-users and irregular 
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Cardiovascular diseases 


Yes 


History of smoking 
Yes 

No 

History of stroke 
Yes 

No 

History of MI 

Yes 


No 


TABLE 4: Hospitalization rates in the three two-group comparisons after balancing and matching 


‘Weltmectin 
non-users (n 


48/4 ,582) (3.0%) 


1/9 (11.1%) 


46/1,533 (3.0%) 


1/6 (16.7%) 


46/1,536 (3.0%) 


0/21 (0.0%) 


4711,521 (3.1%) 


0/2 (0.0%) 


47/1,540 (3.1%) 


0/1 (0.0%) 


47/1,541 (3.0%) 


(dgan 
ivermectin 


RELS h2- 1742) 


1/16 (6.3%) 


37/1,526 (2.4%) 


0/6 (0.0%) 


38/1,536 (2.5%) 


0/23 (0.0%) 


38/1,519 (2.5%) 


0/3 (0.0%) 


38/1,539 (2.5%) 


1/3 (33.3%) 


37/1,539 (2.4%) 


fra Nistad hoa peal risk 
ratio (95% cmand p- 
Wu fp} — 1.27) (0.38) 


0.53 (0.03 — 9.71) (0.67) 


0.80 (0.52 — 1.25) (0.33) 


0.28 (0.01 — 8.42) (0.47) 


0.82 (0.53 — 1.27) (0.38) 


0.91 (0.02 — 48.2) (0.96) 


0.80 (0.52 — 1.24) (0.33) 


0.71 (0.01 — 49.7) (0.88) 


0.80 (0.52 — 1.24) (0.32) 


1.80 (0.04 — 79.4) (0.76) 


0.78 (0.51 — 1.21) (0.27) 


the groups of non-users and irregular ivermectin users 


“Multivariate adjusted hospital 
risk ratio (95% chand p-value 


Oppy (0.49 — 1.12) ( 


n/a (0.99) 


0.70 (0.46 — 1.06) (0.09) 


n/a (0.98) 


0.74 (0.49 — 1.11) (0.14) 


0.97 (n/a) (1.00) 


0.71 (0.48 — 1.07) (0.10) 


n/a (1.00) 


0.72 (0.48 — 1.08) (0.11) 


n/a (0.99) 


0.70 (0.46 — 1.07) (0.09) 


MI = myocardial infarction; COPD = chronic obstructive pulmonary disease; y/o = years old; Cl = confidence interval; n/a = not applicable. 


Propensity score matched 
regular users and irregular 
users 


Overall 

Age 

<30 y/o 

30-50 y/o 

>50 y/o 

Sex 

Female 

Male 

Race 
Afro-Brazilian 
Mixed 
Caucasian 
Asian-Brazilian 
Type 2 diabetes 


Yes 


Regular 
ivermectin 
users (n = 283) 


0/283 (0.0%) 


0/39 (0.0%) 
0/131 (0.0%) 


0/113 (0.0%) 


0/141 (0.0%) 


0/142 (0.0%) 


0/4 (0.0%) 
0/58 (0.0%) 
2/221 (0.0%) 


0/0 


0/9 (0.0%) 
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Irregular 
ivermectin 


users (n = 283) 


10/283 (3.5%) 


0/60 (0.0%) 
3/132 (2.3%) 


7191 (7.7%) 


2/155 (1.3%) 


8/128 (6.3%) 


0/5 (0.0%) 
3/68 (4.4%) 
7/209 (3.3%) 


0/1 (0.0%) 


1/10 (10.0%) 


Unadjusted hospital risk 
ratio (95% Cl) and p- 
value (p) 


0.05 (0.003 — 0.79) (0.034) 


1.53 (0.03 — 78.8) (0.83) 
0.14 (0.01 — 2.75) (0.20) 


0.05 (0.003 — 0.88) (0.041) 


0.22 (0.01 — 4.56) (0.33) 


0.05 (0.003 — 0.87) (0.04) 


1.22 (0.02 = 74.3) (0.92) 
0.16 (0.01 — 3.16) (0.23) 
0.26 (0.05 — 1.28) (0.099) 


3.00 (0.02 — 473.1) (0.67) 


0.33 (0.01 — 9.26) (0.52) 


Multivariate adjusted 
hospital risk ratio (95% Cl) 
and p-value (p) 


0.00 (n/a) (<0.0001) 


1.00 (n/a) (1.00) 
n/a 


n/a (n/a) (<0.0001) 


n/a (n/a) (<0.0001) 


n/a (n/a) (<0.0001) 


1.00 (n/a) (1.00) 
n/a 
n/a (n/a) (<0.0001) 


n/a 
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Cureus 


0/274 (0.09 9/273 (3.3% 0.05 (0.003 — 0.88) (0,0 0.00 <0.0001 
Mo auneny score matched egular 2) EAEE PNA PENEI =) k RETO T 
regular users and irregular ivermectin ivermectin ratio (95% CI) and p- hospital risk ratio (95% CI) 


nsion 


users (n = 283) 


users (n = 283) 


value (p) 


and p-value (p) 


Yes 0/23 (0.0%) 1/20 (5.0%) 0.28 (0.01 — 7.18) (0.44) n/a 

No 0/260 (0.0%) 9/263 (3.4%) 0.05 (0.003 — 0.89) (0.041) n/a (n/a) (<0.0001) 
Asthma 

Yes 0/0 0/0 n/a n/a 

No 0/283 (0.0%) 10/283 (3.5%) 0.05 (0.003 — 0.79) (0.034) 0.00 (n/a) (<0.0001) 
COPD 

Yes 0/1 (0.0%) 0/0 0.33 (0.002 — 52.6) (0.67) n/a 

No 0/282 (0.0%) 10/283 (3.5%) 0.05 (0.003 — 0.79) (0.034) 0.00 (n/a) (<0.0001) 
Other respiratory diseases 

Yes 0/1 (0.0%) 0/0 0.33 (0.002 — 52.6) (0.67) n/a 

No 0/282 (0.0%) 10/283 (3.5%) 0.05 (0.003 — 0.79) (0.034) 0.00 (n/a) (<0.0001) 
Cardiovascular diseases 

Yes 0/2 (0.0%) 0/5 (0.0%) 2.20 (0.03 — 146.0) (0.71) 1.00 (n/a) (1.00) 

No 0/281 (0.0%) 10/278 (3.6%) 0.05 (0.003 — 0.78) (0.033) 0.00 (n/a) (<0.0001) 
Cancer 

Yes 0/2 (0.0%) 0/2 (0.0%) 1.00 (0.01 — 73.3) (1.00) 1.00 (n/a) (1.00) 

No 0/281 (0.0%) 10/281 (3.6%) 0.05 (0.003 — 0.79) (0.034) 0.00 (n/a) (<0.0001) 
History of smoking 

Yes 0/3 (0.0%) 0/1 (0.0%) 0.43 (0.01 — 33.6) (0.70) 1.00 (n/a) (1.00) 

No 0/280 (0.0%) 10/282 (3.5%) 0.05 (0.003 — 0.79) (0.034) 0.00 (n/a) (<0.0001) 
History of stroke 

Yes 0/1 (0.0%) 0/1 (0.0%) 1.00 (0.01 — 92.4) (1.00) n/a 

No 0/282 (0.0%) 10/282 (3.5%) 0.05 (0.003 — 0.79) (0.034) 0.00 (n/a) (<0.0001) 
History of MI 

Yes 0/0 0/0 n/a n/a 

No 0/283 (0.0%) 10/283 (3.5%) 0.05 (0.003 — 0.79) (0.034) 0.00 (n/a) (<0.0001) 


TABLE 5: Hospitalization rates in the three two-group comparisons after balancing and matching 
the groups of regular and irregular ivermectin users 


MI = myocardial infarction; COPD = chronic obstructive pulmonary disease; y/o = years old; CI = confidence interval; n/a = not applicable. 


Figure 3 illustrates differences in hospitalization rates in the overall population between matched groups. 
Balanced and matched groups of non-users and regular users (283 subjects in each group) showed 13 
hospitalizations among non-users (4.6% hospitalization rate) and zero hospitalizations among regular 
users (0.0% hospitalization rate), a 100% reduction after adjustment for variables (RR: 0.00; 95% CI: not 
applicable (n/a); p < 0.0001). Between non-users and irregular users (n = 1,542 in each group), there were 47 
hospitalizations among non-users (3.0% hospitalization rate) and 38 hospitalizations among irregular 
ivermectin users (2.5% hospitalization rate), i.e., a 29% reduction (RR: 0.71; 95% CI: 0.49 - 1.05; p = 0.099). 
Between regular and irregular users (n = 283 in each group), there were 10 hospitalizations among irregular 
users (3.5% hospitalization rate) and zero hospitalizations among regular users (0.0% hospitalization rate), 
i.e., a 100% reduction after adjustment for variables (RR: 0.00; 95% CI: n/a; p < 0.0001). Precise comparisons 
between subpopulations of regular users and non-users and between regular users and irregular users were 
precluded due to a lack of hospitalizations among regular users, as observed in Figure 3. 
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| 
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(RR, 0.00; 95%C, n/a, p < 0.001) (RR, 0.71; 95%CI, 0.49 = 1.05, p = 0.099) (RR, 0.00; 95%0, n/a, p < 0.0001) 


FIGURE 3: Hospitalization rates for the overall population in post- 
matched groups 


RR = risk ratio; Cl = confidence interval. 


Mortality rates among ivermectin non-users, regular users, and 
irregular users 


Supplement Appendix 1B shows mortality rates in ivermectin non-users, regular users, and irregular users 
before matching is described. Tables 6-8 and Figure 4 show mortality rates for each of the three 
combinations of post-matched groups of ivermectin non-users and regular users, non-users and irregular 
users, and regular and irregular users. 


Propensity score Ivermectin Regular Unadjusted mortality risk Multivariate adjusted mortality 
matched non-users and non-users (n ivermectin ratio (95% Cl) and p-value risk ratio (95% Cl) and p-value 
regular users = 283) users (n = 283) (p) (p) 

Overall 15/283 (5.3%) 2/283 (0.7%) 0.13 (0.03 — 0.56) (0.006) 0.08 (0.02 — 0.35) (0.0008) 

Age 

<30 y/o 0/63 (0.0%) 0/39 (0.0%) 1.61 (0.03 — 82.7) (0.81) n/a (1.00) 

30-50 y/o 1/152 (0.7%) 0/131 (0.0%) 0.38 (0.02 — 9.51) (0.56) n/a (0.98) 

>50 y/o 14/68 (20.6%) 2/113 (1.8%) 0.07 (0.02 — 0.32) (0.0006) 0.08 (0.02 — 0.37) (0.001) 

Sex 

Female 8/156 (5.1%) 0/141 (0.0%) 0.06 (0.004 — 1.08) (0.056) 0.00 (0.00 — 0.00) (<0.0001) 

Male 7/127 (5.5%) 2/142 (1.4%) 0.24 (0.05 — 1.20) (0.083) 0.15 (0.03 — 0.70) (0.015) 

Race 

Afro-Brazilian 1/9 (11.1%) 0/4 (0.0%) 0.63 (0.02 — 18.8) (0.79) n/a 

Mixed 1/58 (1.7%) 0/58 (0.0%) 0.33 (0.01 — 8.21) (0.50) n/a 

Caucasian 13/213 (6.1%) 2/221 (0.9%) 0.15 (0.03 — 0.66) (0.012) n/a 

Asian-Brazilian 0/3 (0.0%) 0/0 7.00 (0.05 — 953.3) (0.44) n/a 

Type 2 diabetes 

Yes 3/10 (30.0%) 1/9 (11.1%) 0.29 (0.02 — 3.48) (0.33) 0.33 (0.04 — 2.58) (0.16) 

No 12/273 (4.4%) 1/274 (0.4%) 0.08 (0.01 — 0.62) (0.015) 0.05 (0.01 — 0.37) (0.004) 
Hypertension 

Yes 6/21 (28.5%) 1/23 (4.3%) 0.11 (0.01 — 1.04) (0.054) 0.16 (0.02 — 1.16) (0.07) 

No 9/262 (3.4%) 1/260 (0.4%) 0.11 (0.01 — 0.86) (0.036) 0.06 (0.01 — 0.49) (0.009) 

Asthma 

Yes 0/0 0/0 n/a n/a 
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History of smoking 
Yes 

No 

History of stroke 
Yes 

No 

History of MI 

Yes 


No 


LENGEN 


= 283) 


0/0 


15/283 (5.3%) 


0/1 (0.0%) 


15/282 (5.3%) 


0/1 (0.0%) 


15/282 (5.3%) 


1/2 (50.0%) 


14/281 (5.0%) 


0/2 (0.0%) 


15/281 (5.3%) 


1/1 (100.0%) 


14/282 (5.0%) 


0/0 


15/283 (5.3%) 


Dicis liE] 


users (n = 283) 


0/1 (0.0%) 


2/282 (0.7%) 


0/1 (0.0%) 


2/282 (0.7%) 


0/2 (0.0%) 


2/281 (0.7%) 


0/2 (0.0%) 


2/281 (0.7%) 


0/3 (0.0%) 


2/280 (0.7%) 


0/1 (0.0%) 


2/282 (0.7%) 


0/0 


2/283 (0.7%) 


PSSS Ca hue 


0.33 (0.002 — 52.6) (0.67) 


0.13 (0.03 — 0.56) (0.007) 


1.00 (0.01 — 92.4) (1.00) 


0.13 (0.03 — 0.56) (0.006) 


0.60 (0.007 — 49.5) (0.82) 


0.13 (0.03 — 0.56) (0.007) 


0.20 (0.005 — 8.83) (0.40) 


0.14 (0.03 — 0.61) (0.009) 


0.71 (0.01 — 49.7) (0.88) 


0.13 (0.03 — 0.56) (0.007) 


0.11 (0.001 — 10.3) (0.34) 


0.14 (0.03 — 0.61) (0.009) 


n/a 


0.13 (0.03 — 0.56) (0.007) 


ultivariate adjusted mortality 
fsk ratio (95% Cl) and p-value 
(p) 


n/a 


n/a 


n/a 


n/a 


n/a 


n/a 


n/a 


0.09 (0.02 — 0.37) (0.001) 


n/a (1.00) 


0.08 (0.02 — 0.36) (0.0008) 


0.01 (0.001 — 0.02) (<0.0001) 


0.10 (0.02 — 0.40) (0.001) 


n/a 


TABLE 6: Mortality rates in the three ivermectin two-group matches of non-users and regular 


users 


MI = myocardial infarction; COPD = chronic obstructive pulmonary disease; y/o = years old; Cl = confidence interval; n/a = not applicable. 


Propensity score 
matched non-users and 
irregular users 


Overall 
Age 

<30 y/o 
30-50 y/o 
>50 y/o 
Sex 
Female 
Male 


Race 


Ivermectin 
non-users (n 
= 1,542) 


46/1,542 (3.0%) 


0/410 (0.0%) 
5/808 (0.6%) 


41/324 (12.7%) 


27/846 (3.2%) 


19/696 (2.7%) 
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Irregular 
ivermectin 
users (n = 
1,542) 


29/1,542 (1.9%) 


0/397 (0.0%) 
2/775 (0.3%) 


27/370 (7.3%) 


15/853 (1.8%) 


14/689 (2.0%) 


Unadjusted mortality risk 


ratio (95% Cl) and p- 
value (p) 


0.62 (0.39 — 0.99) (0.049) 


1.03 (0.02 — 52.2) (0.99) 
0.42 (0.08 — 2.15) (0.29) 


0.54 (0.33 — 0.91) (0.019) 


0.54 (0.29 — 1.03) (0.061) 


0.74 (0.37 — 1.49) (0.40) 


Multivariate adjusted 
mortality risk ratio (95% Cl) 
and p-value (p) 


0.63 (0.40 — 0.99) (0.049) 


1.00 (0.63 — 1.59) (1.00) 
0.42 (0.08 — 2.14) (0.30) 


0.58 (0.36 — 0.92) (0.02) 


0.55 (0.30 — 0.99) (0.049) 


0.58 (0.30 — 1.12) (0.11) 
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Cureus 


Afro-Brazilian 


Propensity score 
hed non-users and 
irregular users 


Caucasian 
Asian-Brazilian 
Type 2 diabetes 
Yes 

No 
Hypertension 
Yes 

No 


Asthma 


Other respiratory diseases 
Yes 
No 


Cardiovascular diseases 


Cancer 


No 

History of smoking 
Yes 

No 

History of stroke 
Yes 

No 

History of MI 

Yes 


No 


2/45 (4.4%) 
Ivermectin 


eye 


36/1,114 (3.2%) 


1/32 (3.1%) 


10/37 (27.0%) 


36/1,505 (2.4%) 


16/86 (18.6%) 


30/1,456 (2.1%) 


1/6 (16.7%) 


45/1,536 (2.9%) 


0/1 (0.0%) 


46/1,541 (3.0%) 


1/3 (33.3%) 


45/1,539 (2.9%) 


1/9 (11.1%) 


45/1,533 (2.9%) 


1/6 (16.7%) 


45/1,536 (2.9%) 


1/21 (4.8%) 


45/1,521 (3.0%) 


0/2 (0.0%) 


46/1,540 (3.0%) 


0/1 (0.0%) 


46/1,541 (3.0%) 


0/37 (0.0%) 
rreguiar 


ivermectin 
7hsets (12>) 


1,542) 
21/1,102 (1.9%) 


1/30 (3.3%) 


3/40 (7.5%) 


26/1,502 (1.7%) 


7/96 (7.3%) 


22/1,446 (1.5%) 


1/6 (16.7%) 


28/1,536 (1.8%) 


0/1 (50.0%) 


29/1,541 (1.9%) 


0/3 (0.0%) 


29/1,539 (1.9%) 


0/16 (0.0%) 


29/1,526 (1.9%) 


0/6 (0.0%) 


29/1,536 (1.9%) 


0/23 (0.0%) 


29/1,519 (1.9%) 


0/3 (0.0%) 


29/1,539 (1.9%) 


0/3 (0.0%) 


29/1,539 (1.9%) 


0.23 (0.01 — 4.99) (0.35) 
Unadjusted mortality risk 
į % ȘI 
s SA 
0.58 (0.34 — 1.00) (0.05) 


1.07 (0.06 — 17.9) (0.96) 


0.22 (0.05 — 0.87) (0.031) 


0.72 (0.43 — 1.20) (0.20) 


0.34 (0.13 — 0.88) (0.026) 


0.73 (0.42 — 1.27) (0.26) 


1.00 (0.05 — 20.8) (1.00) 


0.62 (0.38 — 0.99) (0.046) 


1.00 (0.01 — 92.4) (1.00) 


0.62 (0.39 — 0.99) (0.049) 


0.24 (0.01 — 8.62) (0.43) 


0.64 (0.40 — 1.02) (0.062) 


0.17 (0.01 — 4.68) (0.30) 


0.64 (0.40 — 1.03) (0.064) 


0.28 (0.01 — 8.42) (0.47) 


0.64 (0.40 — 1.02) (0.062) 


0.29 (0.01 — 7.54) (0.46) 


0.64 (0.40 — 1.02) (0.063) 


0.71 (0.01 — 49.7) (0.88) 


0.62 (0.39 — 0.99) (0.049) 


0.43 (0.01 — 33.6) (0.70) 


0.62 (0.39 — 0.99) (0.049) 


n/a 
Multivariate adjusted 


lity ri ig {95% CI 
pertela riae. C 
0.53 (0.32 — 0.89) (0.016) 


0.81 (0.07 — 9.99) (0.87) 


0.32 (0.10 — 1.04) (0.057) 


0.64 (0.39 — 1.04) (0.069) 


0.38 (0.17 — 0.87) (0.022) 


0.66 (0.39 — 1.12) (0.12) 


3.95 (0.02 — 789.8) (0.61) 


0.56 (0.36 — 0.88) (0.011) 


1.00 (0.64 — 1.56) (1.00) 


0.56 (0.36 — 0.88) (0.011) 


0.06 (0 — 4178.4) (0.62) 


0.58 (0.37 — 0.91) (0.017) 


0.04 (0.03 — 0.07) (<0.0001) 


0.57 (0.36 — 0.88) (0.012) 


n/a 


0.58 (0.37 — 0.90) (0.016) 


0.00 (<0.0001) 


0.57 (0.37 — 0.90) (0.015) 


0.40 (0.26 — 0.62) (<0.0001) 


0.56 (0.36 — 0.88) (0.011) 


0.04 (0.03 — 0.07) (<0.0001) 


0.57 (0.36 — 0.88) (0.012) 


TABLE 7: Mortality rates in the three ivermectin two-group matches of non-users and irregular 


users 


MI = myocardial infarction; COPD = chronic obstructive pulmonary disease; y/o = years old; Cl = confidence interval; n/a = not applicable. 
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Cureus 


Propensity score matched 
regular users and irregular 


users 


Overall 

Age 

<30 y/o 

30-50 y/o 

>50 y/o 

Sex 

Female 

Male 

Race 
Afro-Brazilian 
Mixed 
Caucasian 
Asian-Brazilian 
Type 2 diabetes 
Yes 

No 
Hypertension 


Yes 


Other respiratory diseases 
Yes 
No 


Cardiovascular diseases 


Cancer 


No 
History of smoking 
Yes 


No 


Regular 
ivermectin 
users (n = 283) 


2/283 (0.7%) 


0/39 (0.0%) 
0/131 (0.0%) 


2/113 (1.8%) 


0/141 (0.0%) 


2/142 (1.4%) 


0/4 (0.0%) 
0/58 (0.0%) 
2/221 (0.9%) 


0/0 


1/9 (11.1%) 


1/274 (0.4%) 


1/23 (4.3%) 


1/260 (0.4%) 


0/0 


2/283 (0.7%) 


0/1 (0.0%) 


2/282 (0.7%) 


0/1 (0.0%) 


2/282 (0.7%) 


0/2 (0.0%) 


2/281 (0.7%) 


0/2 (0.0%) 


2/281 (0.7%) 


0/3 (0.0%) 


2/280 (0.7%) 
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Irregular 
ivermectin 
users (n = 283) 


10/283 (3.5%) 


0/60 (0.0%) 
0/132 (0.0%) 


10/91 (11.0%) 


4/155 (2.6%) 


6/128 (4.7%) 


0/5 (0.0%) 
3/68 (4.4%) 
7/209 (3.3%) 


0/1 (0.0%) 


1/10 (10.0%) 


9/273 (3.3%) 


1/20 (5.0%) 


9/263 (3.4%) 


0/0 


10/283 (3.5%) 


0/0 


10/283 (3.5%) 


0/0 


10/283 (3.5%) 


0/5 (0.0%) 


10/278 (3.6%) 


0/2 (0.0%) 


10/281 (3.6%) 


0/1 (0.0%) 


10/282 (3.5%) 


Unadjusted mortality 


risk ratio (95% Cl) and p- 


value (p) 


0.19 (0.04 — 0.89) (0.036) 


1.53 (0.03 - 78.8) (0.83) 
1.01 (0.02 — 51.2) (1.00) 


0.15 (0.03 — 0.68) (0.015) 


0.12 (0.01 — 2.23) (0.15) 


0.29 (0.06 — 1.47) (0.13) 


1.22 (0.02 — 74.7) (0.92) 
0.16 (0.01 — 3.16) (0.23) 
0.26 (0.05 — 1.28) (0.099) 


3.00 (0.02 — 473.1) (0.67) 


1.13 (0.06 — 21.1) (0.94) 


0.11 (0.01 — 0.85) (0.035) 


0.86 (0.05 — 14.8) (0.92) 


0.11 (0.01 — 0.79) (0.036) 


n/a 


0.19 (0.04 — 0.89) (0.036) 


0.33 (0.002 — 52.6) (0.67) 


0.19 (0.04 — 0.90) (0.036) 


n/a 


0.19 (0.04 — 0.90) (0.036) 


2.20 (0.03 — 146.1) (0.71) 


0.19 (0.04 — 0.88) (0.034) 


1.00 (0.01 — 73.3) (1.00) 


0.19 (0.04 — 0.89) (0.036) 


0.43 (0.01 — 33.6) (0.70) 


0.20 (0.04 — 0.90) (0.036) 


Multivariate adjusted 


mortality risk ratio (95% Cl) 


and p-value (p) 


0.16 (0.04 — 0.71) (0.016) 


1.00 (0.22 — 4.46) (1.00) 
1.00 (0.22 — 4.46) (1.00) 


0.16 (0.04 — 0.72) (0.017) 


0.00 (n/a) (0.98) 


0.25 (0.05 — 1.19) (0.082) 


n/a 
n/a 
n/a 


n/a 


0.88 (0.07 — 11.6) (0.92) 


0.09 (0.01 — 0.69) (0.021) 


0.94 (0.06 — 13.9) (0.96) 


0.09 (0.01 — 0.67) (0.019) 


n/a 


0.16 (0.04 — 0.71) (0.016) 


n/a 


0.16 (0.04 — 0.72) (0.017) 


n/a 


0.16 (0.04 — 0.72) (0.017) 


0.52 (0.11 — 2.30) (0.38) 


0.16 (0.04 — 0.71) (0.016) 


1.00 (0.22 — 4.46) (1.00) 


0.16 (0.04 — 0.70) (0.016) 


n/a 


0.16 (0.03 — 0.72) (0.017) 
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Cureus 


HBSS AEA ore matched Regular Irregular Unadjusted mortality Multivariate adjusted 
regular users and irregular ivermectin ivermectin risk ratio (95% Cl) and p- mortality risk ratio (95% Cl) 

Yusers OvseODa)= 283) OvBEDDfO)= 283) 1vatuo (p) — 92.4) (1.00) 100 (P-2alue ($8) (1.00) 

No 2/282 (0.7%) 10/282 (3.5%) 0.19 (0.04 — 0.89) (0.036) 0.16 (0.04 — 0.72) (0.017) 

History of MI 

Yes 0/0 0/0 n/a n/a 

No 2/283 (0.7%) 10/283 (3.5%) 0.19 (0.04 — 0.89) (0.036) 0.16 (0.04 — 0.71) (0.016) 


TABLE 8: Mortality rates in the three ivermectin two-group matches of regular users and irregular 
users 


MI = myocardial infarction; COPD = chronic obstructive pulmonary disease; y/o = years old; CI = confidence interval; n/a = not applicable. 
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PROPENSITY SCORE MATCHING 
MATCHED 
NON-USE AND REGULAR USE NON-USE AND IRREGULAR USE REGULAR AND IRREGULAR USE GROUPS 
NON IVERMECTIN REGULAR NON IVERMECTIN IRREGULAR REGULAR IRREGULAR 
USE IVERMECTIN USE USE IVERMECTIN USE IVERMECTIN USE IVERMECTIN USE COVID-19 
l = mortali 
283 subjects 283 subjects | 1,452 subjects 1,542 subjects 283 subjects 283 subjects rate ty 
i | | ! | 
| 15 deaths 02 deaths | 46 deaths 29 deaths 02 deaths 10 deaths 
5.3% 0.7% 3.0% 1.9% 0.7% 3.5% 
mortality rate mortality rate mortality rate mortality rate mortality rate mortality rate Overall 
J L population 
92% vs non-users 37% vs non-users 84% vs irregular users 
(RR, 0.08; 95%C1, 0.02 - 0.35; p = 0.0008) (RR, 0.63; 95%CI, 0.40 - 0.99; p = 0.049) IRR, 0.16 95%CI, 0.04 — 0.71; p = 0.016) 
08 deaths O deaths 27 deaths ] 15 deaths O deaths 04 deaths 
(out of 156) (out of 141) (out of 846) (out of 853) (out of 141) (out of 155) 
5.1% 0.0% 3.2% 1.8% 0.0% 2.6% 
mortality rate mortality rate mortality rate mortality rate mortality rate mortality rate) Females 
100% vs non-users 145% vs non-users 100% vs irregular users 
(RR, 0.00; 95%C1, n/a; p < 0.0001) (RR, 0.55; 95%, 0.30 - 0.99; p = 0.049) (RR, 0.00; 95%CI, n/a; p = 0.98) 
O7deaths 02 deaths 19deaths | 14 deaths 02 deaths 06 deaths 
(out of 127) {out of 142) (out of 696) | (out of 689) (out of 142) (out of 128) 
5.5% 1.4% 2.7% 2.0% 1.4% 4.7% 
mortality rate mortality rate mortality rate mortality rate mortality rate mortality rate M al es 
485% vs non-users 442% vs non-users 75% vs irregular users 
(RR, 0.15; 95%, 0.03 — 0.70; p = 0.015) (RR, 0.58; 95%CI, 0.30 — 1.12; p = 0.11) (RR, 0.25 95%), 0.05 ~ 1.19; p = 0.082) 
27 death: 02 deaths 10 deaths 
en E) na EREA (out of 5 (out of 113) {out of 91) 
20.6% 1.8% 12.7% 7.3% 1.8% 11.0% 
mortality rate mortality rate mortality rate mortality rate mortality rate mortality rate Above 50 y/o 
92% vs non-users 42% vs non-users 84% vs irregular users 
(RR, 0.08; 95%), 0.02 0.37; p = 0.001) (RR, 0.58; 95%C1, 0.36 - 0.92; p = 0.02) (RR, 0.16 95%C1, 0.04 — 0.72; p = 0.017) 
] O1death 03 death: O1death Oi death 
[ Gem I eee RN (out of 40) (out of 03) (out of 10) 
30.0% 11.1% 27.0% 7.5% 11.1% 10.0% Type 2 
mortality rate mortality rate mortality rate mortality rate mortality rate mortality rate diabetes 
67% vs non-users 68% vs non-users 112% vs irregular users (adjusted) 
(RR, 0.33; 95%CI, 0.04 ~ 2.58; p = 0.16) (RR, 0.32; 95%C1, 0.10 — 1.04; p = 0.057) (RR, 0.88; 95%CI, 0.07 — 11.6; p = 0.92) 
| 06 deaths 16 deaths 07 deaths Old 
| {out of 21) (out of 86) (out of 96) (out of 23) 
28.5% 4.3% 18.6% 7.3% 5.0% 
mortality rate mortality rate mortality rate mortality rate mortality rate mortality rate Hypertension 
84% vs non-users 62% vs non-users 6% vs irregular users 
(RR, 0.16; 95%0, 0.02 — 1.16; p = 0.07) (RR, 0.38; 95%CI, 0.17 -0.87; p = 0.022) (RR, 0.94; 95%C1, 0.06 - 13.9; p = 0.96) 


FIGURE 4: Mortality rates in the post-matched, overall population, and 
subpopulation groups 


RR = risk ratio; CI = confidence interval. 


Between matched groups of non-users and regular users (n = 283 in each group), mortality rate was 5.3% (15 
deaths) among non-users and 0.7% (two deaths) among regular users; a 92% reduction in mortality rate (RR: 
0.08; 95% CI: 0.02 - 0.35; p = 0.00083). Compared to non-users, reductions in mortality rate among regular 
users were 100% among females (eight deaths among 156 non-users and zero deaths among 141 regular 
users; RR: 0.00; 95% CI: n/a; p < 0.0001) and 85% among males (seven deaths among 127 non-users and two 
deaths among 142 regular users; RR: 0.15; 95% CI: 0.03 - 0.70; p = 0.015) with 92% for subjects above 50 
years of age (14 deaths among 68 non-users and two deaths among 113 regular users; RR: 0.08; 95% CI: 0.02 
- 0.37; p = 0.001). There was statistically a non-significant 67% reduction for T2D (three deaths among 10 
non-users and one death among nine regular users; RR: 0.33; 95% CI: 0.04 - 2.58; p = 0.16), and 84% among 
subjects with hypertension (six deaths among 21 non-users and one death among 23 regular users; RR: 0.16; 
95% CI: 0.02 - 1.16; p = 0.07). 


Between matched groups of non-users and irregular users (n = 1,542 in each group), there was a 3.0% 
mortality rate (46 deaths) among non-users and a 1.9% mortality rate (29 deaths) among irregular users, 
showing a 37% reduction in mortality rate (RR compared to non-users: 0.63; 95% CI: 0.40 - 0.99; p = 0.049). 
A 45% reduction in mortality rate occurred among females; 3.2% for non-users (27 death among 846) and 
1.8% for irregular users (15 deaths among 853) (RR: 0.55; 95% CI: 0.30 - 0.99; p = 0.049); and 42% reduction 
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Cureus 


occurred for males; 2.7% of non-users (19 deaths among 696) and 2.0% of irregular users (14 deaths among 
689) (RR: 0.58; 95% CI: 0.30 - 1.12; p = 0.11). Mortality rate for subjects above 50 years old was 12.7% for 324 
non-users (41 deaths) and 7.3% for 370 regular users (27 deaths); a 42% reduction in mortality rate (RR: 
0.58; 95% CI: 0.36 - 0.92; p = 0.02). Participants with T2D had 27.0% mortality rate for 37 non-users (10 
deaths) and 7.5% for 40 irregular users (three deaths); a 68% reduction in mortality rate among participants 
with T2D (RR: 0.32; 95% CI: 0.10 - 1.04; p = 0.057). Those with hypertension had a 62% reduction in 
mortality rate; 18.6% of 86 non-users (16 deaths) and 7.3% of 96 irregular users (seven deaths) (RR: 0.38; 
95% CI: 0.17 - 0.87; p = 0.022). In sub-populations without comorbidities, reductions in mortality rates were 
between 40% and 45%. 


When groups of regular users and irregular users are matched (283 subjects in each group), there was a 0.7% 
and 3.5% (two deaths and 10 deaths) mortality rate among regular and irregular users, reflecting a 
multivariate-adjusted 84% reduction in mortality rate (RR: 0.16; 95% CI: 0.04 - 0.71; p = 0.016). The small 
number of events between these two groups precludes more statistically significant differences, despite 
large effect size and differences, in particular in subgroups with fewer subjects. The mortality rate was 2.6% 
(four deaths out of 155) among non-user females and 0.0% (out of 141 females) among regular user females. 
There was a 4.7% mortality rate (six deaths) among 128 non-user males and a 1.4% mortality rate (two 
deaths) among 142 regular user males, showing a reduction of 75% (RR: 0.25; 95% CI: 0.05 - 1.19; p = 0.082) 
in mortality rate. Reduction in mortality rate was 84% for those over 50 years of age; 11.0% for non-users 
(10 deaths among 91) and 1.8% for regular users (two deaths among 113) (RR: 0.16; 95% CI: 0.04 - 0.72; p= 
0.017). Among participants with T2D, mortality rate was 10.0% (one death) among 10 irregular users and 
11.1% (one death) among nine regular users, which is statistically similar between groups (RR: 0.88; 95% CI: 
0.07 - 11.6; p = 0.92). Subjects with hypertension had 5.0% mortality rate (one death) among 20 irregular 
users and 4.3% mortality rate (one death) among 23 regular users, which is similar between groups (RR: 
0.94; 95% CI: 0.06 - 13.9; p = 0.96). 


Risk of dying from COVID-19 between ivermectin non-users, regular 
users, and irregular users 


Considering the population and participants of Itajai, as well as inhabitants of Itajai, who did not use 
ivermectin prophylactically, the unadjusted risk of dying from COVID-19 was 1,730 in every 1,000,000 
subjects among non-users, 240 among regular users, and 850 among irregular users. Compared to non-users, 
the risk of dying from COVID-19 was 86% lower in regular users (RR: 0.14; 95% CI: 0.03 - 0.57; p = 0.006) and 
51% lower in irregular users (RR: 0.49; 95% CI: 0.32 - 0.76; p = 0.001). The risk of dying from COVID-19 was 
72% lower in regular users than irregular users (RR: 0.28; 95% CI: 0.07 - 1.18; p = 0.089). Figure 5 illustrates 
the risk of dying from COVID-19 in each population. 


2022 Kerr et al. Cureus 14(8): e28624. DOI 10.7759/cureus.28624 20 of 37 


Cureus 


Citywide program of COVID-19 prophylaxis with ivermectin (Itajaí, Brazil) - 
Risk of dying from COVID-19 


NON IVERMECTIN USE MEES MEE 
45,716 8,325 33,971 
OVERALL 


79 deaths 2 deaths 29 deaths 


sh 


0.0173% 0.00024% 0.0085% 
risk of death = risk of death = risk of death = 
1,730/1,000,000 240/1,000,000 850/1,000,000 


86% reduction 51% reduction 
vs non-use vs non-use 
0.14 (0.03 — 0.57) 0.49 (0.32 — 0.76) 
p =0.006 p=0.001 


72% reduction 
vs irregular use 
0.28 (0.07 — 1.18) 
p = 0.083 


FIGURE 5: Risk of dying from COVID-19 among ivermectin non-users, 
regular users, and irregular users 


RR = risk ratio; CI = confidence interval. 


Discussion 
The program in Itajaí, Brazil: ivermectin prophylaxis for COVID-19 


The present study provides in-depth results on the prospective study of ivermectin as prophylaxis for 
COVID-19 in Itajaí, located in Southern Brazil. Particularities of Itajaí included its dynamic population due 
to the presence of an overwhelmingly large port compared to the size of the city. This explained why the city 
was one of the first in the state to reach 1,000 cases in 2020 [26]. In the past, the city experienced some of 
the highest rates of HIV infections in Brazil [27], partially substantiated by being a port city, an 
"independent" predictor of a higher prevalence of HIV infection [28]. 


The decision to adopt a prophylaxis program with ivermectin in Itajaí was based on (1) the fact that case 
numbers rose rapidly and at a higher speed than in other cities; (2) the inability to isolate port workers in the 
absence of pharmacological or non-pharmacological therapies for COVID-19; (3) because it had already 
been proven to be a potent antiviral for over 20 viruses, studied independently and peer-reviewed, including 
the first severe acute respiratory syndrome coronavirus (SARS-CoV) epidemic before the COVID-19 
pandemic; and (4) the extensive safety profile and favorable cost-effectiveness of ivermectin. Hence, the 
program of Itajaí strictly followed all bioethical principles using ivermectin as prophylaxis for COVID-19. 
The ivermectin was offered optionally, as prophylaxis for COVID-19, following medical screening by medical 
doctors. 


Ivermectin as a defense against all major COVID-19 outcomes: does it 
depend on the regularity of ivermectin use? 


In our first paper [25], ivermectin was shown to be associated with significant reductions in infection rate 
(44%), hospitalization rate (56%), and mortality rate (68%), when compared to subjects that did not use 
ivermectin prophylactically and irrespective of the regularity of ivermectin use. 
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This study paper analyzes the impact of the regular use of ivermectin on COVID-19 infection. This impact 
included non-users, regular users, and irregular ivermectin users. These groups were estimated from the 
matched population in the city of Itajaí, with an impressive 100% of the population of Itajaí being 
digitalized in the government data system. Their COVID-19 cases, hospitalizations in public hospitals, and 
all deaths due to COVID-19 were strictly followed and recorded. Figure 6 summarizes an overall view of the 


findings of this study. 
NON-PARTICIPANTS PARTICIPANTS 
(n = 45,716) (n = 113,844) 
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een USE] REGULAR USE | (INTERMEDIATE), | IRREGULAR USE | 
| 45,716 8,325 71,548 33,971 
l [ =/> 180mg | ae nating 
in analysis 
45,716 — 8325 a 
2 / 283 N fe 4 542 NS 
/ cases / \ 
25 3,034 / oy \ 
axe Cad | | 3.40% \ | 4.54% \ 
3 Md g f \ irecte | | infection rate | 
oO \ 6.64% J j \ / 
£ infection rate ko 149% JA \ {32% / 
.— | (All population) Sle Soe 
g _ 
= á 3034 O 8 1,842 
ro) [e] d N ta (0) A 38 > 
5 f \ Ver A \ 
= o 0 / 99 \ hospitalizations / hospitalizations \ 
r i \ \ \ 
S| sas = 0.7% | 25% | 
Oo 3 S e| 3.3% \ hospitalization rate j | hospitalization rate | 
" gE \ : / \ J 
a \ hospitalization rate  / \ / 
wj) e N ” sa \ J32% 
ie} N d v8 ee” È ALAL á 
< £ — == 
2 (Among OVID-19 
È infections) 3,034 as a 1,542 
a l A OTN aoe. 
D 79 \ / deaths \ / deaths \ 
2 ami | | 07% \ [| 19% \ 
S Cs 2.6% | mortality rate | | mortality rate | 
Q9 o iny te \ / \ H 
° £ mortality ra \ 185% 7 136% / 
(Among OVID-19 N Ys nonusers“ N E 
infections) — ~ m ——— 
33,971 
aan <Ja 8,325 — 
d ue" y wn yp 2 N 
/ N / deaths À / deaths N 
2 / 79 - | aN 
5 Da | deaths | 0.0024% | | 0.0085% | 
~ 58 | 0 01 73% risk of death J \ risk of death | 
x O 7 H o } \ J 
ka E \ risk of death / ee \ 15t% A 
ÈE Ț 4 Ys romsers Aaa 
(All population) Se i “Sa __ 


All analysis were before matching groups. Analyses of COVID-19 infection rate and risk of dying from COVID-1 
were not adjusted. Analyses of COVID-19 hospitalization and mortality rates were based multivariate 


FIGURE 6: COVID-19 infection, hospitalization, mortality rates, and risk 
of dying from COVID-19 across different patterns of ivermectin use 


This reduction of COVID-19 infection had a significant effect on the reduction of transmission and 
perpetuation of the pandemic in Itajai. Also, the reduction of related hospitalizations and mortality is 
indisputably meaningful. They reduced not only costs and pressure on the health system but saved many 


lives. 


Ivermectin regular users were older (average age = 47 years) compared to irregular users (average age = 41 
years). The non-users (average age = 39.8 years) had approximately 20% to 50% higher prevalence of T2D 
and hypertension. If ivermectin did not work, one would expect higher hospitalizations and mortality rates 
in the group of regular users, which did not happen, as seen in the pre-matched analysis, in Supplement 


Appendix 1. 


Notably, there were no hospitalizations for any of the 289 regular users. After observing matching between 
groups, the reduction in hospitalization rate was 100% in regular users compared to non-users and irregular 
users. Analysis of sub-populations in these two comparisons was unfeasible due to the lack of 
hospitalizations for the regular users. Statistically significant reductions were observed in the 
hospitalization rate for irregular users when compared to non-users (35% reduction; p = 0.03), which was 
more relevant in high-risk populations. This included subjects 50 years of age and above (reduction of 38%; 
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p = 0.027) and those with comorbidities. A 69% reduction was seen among subjects with T2D (p = 0.063), 45% 
among subjects with hypertension (p = 0.10), and 73% among subjects with cardiovascular diseases (p = 
0.23), with reductions similar between males and females. This means that even with uncontrolled, irregular 
use of ivermectin, there is a significant reduction in the number of hospitalizations in COVID-19-infected 
participants. 


The regularity of ivermectin intake demonstrated a progressive impact on the reduction of mortality rate, 
which was more clearly observed after matching groups. The regular users showed a 90% mortality rate 
reduction compared to non-users (p = 0.003) and 79% reduction compared to irregular users (p = 0.05). 
Irregular users had a reduction of 37% compared to non-users (p = 0.63). Reductions among regular users 
were similar (between 86% and 89%) across different high-risk populations (50 years old and above with 
comorbidities). High-risk populations of irregular users had reductions in mortality rate between 34% and 
60% compared to non-users. The most profoundly significant results were for women who used ivermectin 
regularly, with no deaths among all 144 participants. 


The risk of dying from COVID-19, when considering the whole population, was notably lower among regular 
users, compared to both non-users (86% reduction) and irregular users (72% reduction). This risk was also 
lower among irregular users compared to non-users (51% reduction). Since baseline characteristics were not 
present for non-user, non-infected subjects, there were no adjustments to be done for variables relative to 
their chances of dying from COVID-19. 


In common, all outcomes related to COVID-19 infection demonstrated a dose-related response effect, with 
greater reductions in all outcomes with the higher ivermectin intake. This strong correlation reinforces the 
causal relationship between ivermectin intake and protection from COVID-19. Also, although regular users 
still had COVID-19 cases (with a lower infection rate than non-users), these cases tended to be milder, 
compared to non-users or irregular users, as observed in the significant absence of hospitalizations and 
deaths. 


Mechanistically, the accumulated dose of ivermectin, consequently obtained with the regular use of 
ivermectin, had a strong impact on COVID-19-related outcomes, i.e., once infected, higher amounts of 
ivermectin administered related to a better prognosis. Of note, the strict control of which days ivermectin 
was used did not affect the results. 


Although a demonstrative dose-response was observed consistently across the groups (non-users, regular 
users, and irregular users) unexpectedly, the risk of COVID-19 infection was not largely influenced by the 
regularity of ivermectin use (Figure 2). The possible long-term actions of ivermectin, that go beyond its 
serum or cytoplasmatic concentration, may explain the progressive protection with a higher regularity of 
ivermectin use. 


Our results demonstrated protection against COVID-19 when regularly used for two days, every 15 days 
regimen. This prophylactic treatment regimen respected the already extensively known safety profile of 
ivermectin, since, notably, it did not surpass the usual doses for scabies. 


Noteworthy aspects of the study 


Regularity is defined as something happening repeatedly in a fixed pattern. As such, this study determined 
the criteria for regularity to be more than 30 tablets of ivermectin over five months, with a continuous 
supply of ivermectin, determined by the number of tablets prescribed and taken every other week over 12 
weeks. 


To determine different outcomes, it was critical that a correct baseline population was established for each 
outcome. Because there were more than 8,000 subjects from outside the city of Itajaí that participated in the 
study, the infection rate could not be calculated based on the participating subjects because COVID-19 cases 
from other cities were underreported in Itajai among ivermectin non-users. In fact, the “infection rate” of 
overall participants (1.40%) among subjects from other cities (177 cases out of 8,352 subjects) was much 
lower than the infection rates within the city of Itajai. This clearly demonstrated underreporting. 
Calculations were based on participants from Itajaí only, for which COVID-19 cases were strictly controlled. 
Correspondingly, the risk of dying from COVID-19 aims to evaluate the risk of an undesired outcome 
irrespective of how many cases occurred, unlike the mortality rate that included the full population. 


The use of ivermectin was able to reduce COVID-19 infection significantly. A small portion of regular users 
was sufficient to positively affect the city’s numbers related to COVID-19. Unfortunately, because most of 
the population failed to continue in a program of prophylactic ivermectin use, the rise in cases after July 7, 
2020, in the state of Santa Catarina, led to a skewed perception to potentially discredit the efficacy of 
ivermectin. However, misleading this perception, a committed program of ivermectin could have led to a 
huge positive health impact across the whole state. 


Unexpectedly, the different regularity of ivermectin use did not show significant changes in the reduction of 
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COVID-19 infections. One could speculate that subjects that did not obtain ivermectin from the program in 
a regular manner may have acquired ivermectin over the counter, where it was available. However, during 
the first two months of the program, Brazil experienced not only a temporary shortage of ivermectin due to a 
sudden increase in demand, but required a medical prescription and experienced an associated price 
increase by five times, precluding its use outside the program. More importantly, while infection rates did 
not reduce with regular use of ivermectin, compared to irregular users, hospitalization and mortality rates 
reduced substantially, showing a dose-effect response of ivermectin for COVID-19-related outcomes. 


The apparent contradictory lack of hospitalizations, while there were two deaths in the group of regular 
users, may be explained by the fact that patients either used a private hospital outside the city of Itajaí or in 
an institution that was not a hospital. Death reports are mandatorily for public and private hospitals; 
however, hospitalizations are not reported. Another hypothesis is that these deaths occurred without 
hospitalization. Depending on the characteristics and social context of these participants, this is not unusual 
when hospitals get overwhelmed, or when patients avoid seeking hospital care for other reasons [29]. 


Limitations 


Updated medical histories were done for ivermectin users at follow-up appointments with medical doctors 
from the SUS. Regarding the non-users, the participants did not have follow-ups to update their medical 
records. Depending on the calculation methods performed for infection rates, this could create some 
differences. Imprecisions and modifications evident, although minimal, between the first manuscript [25] 
and this study, did not impact the fact that ivermectin use reduced COVID-19-related outcomes. In addition, 
in the present analysis, we did not control for the COVID-19 infection dates. Of note, although there were 
no other hospitals in Itajai, due to the limited capacity of the city hospital, some patients with health 
insurance were transferred to private hospitals outside of Itajai, while some patients without private 
insurance were cared for in institutions that were not hospitals. Unlike hospitalizations, deaths were 
mandatorily reported, which precluded any imprecision in the calculations of the mortality rate. 


The number of tablets was calculated according to body weight. Most of the population used between two 
and three tablets daily for two days, every 15 days. Due to the minimal difference between the number of 
ivermectin tablets used, the amount used (frequency of its use) could be determined with a reasonable level 
of precision. 


This observational study obtained results that presented a high level of certainty by employing strict control 
of the data outcome among COVID-19 cases and strict control of the number of deaths due to COVID-19 in 
the overall population. The fact that PSM was employed for outcomes in such a large population makes 
these data reliable, being sourced from official government databases (datasets: https://osf.io/uxhaf/). 


Final discussion 


Regular use of ivermectin led to a 100% reduction in hospitalization rate, a 92% reduction in mortality rate, 
and an 86% reduction in the risk of dying from COVID-19 when compared to non-users. Irregular use of 
ivermectin led to a 51% reduction in the risk of dying, a 29% reduction in hospitalization rate, and a 37% 
reduction in mortality rate from COVID-19. Statistically significant reductions in hospitalization (100%) and 
mortality rates (84%), and risk of dying from COVID-19 (72%) were observed in regular users when 
compared to irregular users. The response pattern of ivermectin use and level of protection from COVID-19- 
related outcomes was identified and consistent across dose-related levels. The reduction in COVID-19 
infection rate occurred in a consistent and significant dose-dependent manner, with reductions of 49% and 
32% in regular and irregular users, when compared to non-users. The most striking evidence of ivermectin's 
effectiveness was the 100% reduction in mortality for female regular users. 


The analysis of the data gathered from official government databases showed that ivermectin had an 
impactful reduction in the incidence of COVID-19 infection, in a dose-response manner. Even for irregular 
users, benefits were observed. 


The data conclusively show that the risk of dying from COVID-19 was lower for all regular and irregular 
users of ivermectin, compared to non-users, considering the whole population. 


A progressive, dose-response pattern of protection from COVID-19-related outcomes was observed and 
consistent across all levels of ivermectin used. Consequently, the findings in this study show how the risk of 
contracting COVID-19 infection was not greatly influenced by the regularity (regular user = 3.4%; irregular 
user = 4.54%) of ivermectin use, making it very significant as a preventive therapy for COVID-19. 


Finally, the evidence, in this study, added to the efficacy of ivermectin as prophylaxis for COVID-19. There 
are no equivalent RCTs when it comes to the effects of prophylaxis since this was an observational study of a 
strictly controlled population with a great level of control for confounding factors at a magnitude unfeasible 
to be conducted in an RCT. This study demonstrated the effects of ivermectin in real life in an 
overwhelmingly precise manner, close to post-RCT real-life studies [30-32]. The evidence provided by the 
present study is among the strongest and most conclusive data regarding ivermectin efficacy. 
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Female 
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Conclusions 


The regular use of ivermectin decreased hospitalization for COVID-19 by 100%, mortality by 92%, and the 
risk of dying from COVID-19 by 86% when compared to non-users. 


Protection from COVID-19-related outcomes was observed across all levels of ivermectin use, with a notable 
reduction in risk of death in the over 50-year-old population and those with comorbidities. The reduction in 
infection rate was significant, irrespective of the level of ivermectin use. The results of this prospective 
observational study of a strictly controlled population of 223,128 participants reinforce the efficacy of 


ivermectin and the demonstration of a dose-response effect. 


Appendices 
Supplement appendix 1 


Regular use of ivermectin as prophylaxis for COVID-19 led up to a 92% reduction in COVID-19 mortality 
rate in a dose-response manner: results of a prospective observational study of a strictly controlled 


population of 88,012 subjects. 


Index 


Appendix 1A: Hospitalization rates in non-users, irregular users, and regular users before matching 


Appendix 1B: Mortality rates in non-users, irregular users, and regular users before matching 


Appendix 1A: hospitalization rates in non-users, irregular users, and 


regular users before matching 


Tables 9, 10 describe hospitalization rates in ivermectin non-users, irregular users, and strictly regular users, 
before matching groups, and head-to-head comparisons between groups, and Figure 7 illustrates the 
findings in overall and relevant subpopulations. Before matching, there were no hospitalizations among 
ivermectin regular users, 38 hospitalizations among irregular users (2.5%), and 99 hospitalizations among 
non-users (3.3%). After adjustment for variables, reduction in hospitalization rate in ivermectin regular 
users was 98% compared to non-users (RR: 0.02; 95% CI: 0.00 - 0.27; p < 0.0001) and 97% compared to 
irregular users (RR: 0.03; 95% CI: 0.00 - 0.90; p < 0.043). 


Non-matched (pre-matching) groups 


Non-users - 
Overall (n ee 

hospitalizations 
= 4,859) 

(%) (n = 3,034) 
137/4,859 

99/3,034 (3.3%) 
(2.8%) 
2/1,280 

1/844 (0.1%) 
(0.2%) 
12/2,488 

24/1582 (1.5%) 
(1.2%) 
104/1,091 

74/608 (12.7%) 
(9.5%) 
59/2,618 

42/1,624 (2.6%) 
(2.3%) 
68/2,241 

57/1,410 (4.0%) 
(3.0%) 
4/141 

4/100 (4.0%) 
(2.8%) 


Regular users - 
hospitalizations 
(%) (n = 283) 


0/283 (0.0%) 


0/39 (0.0%) 


0/131 (0.0%) 


0/113 (0.0%) 


0/141 (0.0%) 


0/142 (0.0%) 


0/4 (0.0%) 


Irregular users - 
hospitalizations 
(%) (n = 1,542) 


38/1,542 (2.5%) 


1/397 (0.3%) 


71775 (0.9%) 


30/370 (8.1%) 


17/853 (2.0%) 


21/689 (3.0%) 


0/37 (0.0%) 


Overall p-values 


Unadjusted 
overall p- 
value 


<0.0001 


0.46 


<0.0001 


<0.0001 


0.004 


<0.0001 


n/a 


Unadjusted overall p-value (p) and risk ratio 
(95% Cl) between subpopulations 


<0.0001 


3.10 (0.69 — 13.9) (0.14) vs <30 y/o 


67.3 (16.6 — 273.5) (<0.0001) vs <30 y/o / 21.7 
(11.9 - 39.7) (<0.0001) vs 30-50 y/o 


1.36 (0.95 — 1.93) (0.09) vs female 


0.79 


0.98 (0.35 — 2.69) (0.96) vs Caucasians / 1.05 
(0.37 — 3.04) (0.92) vs mixed / 2.60 (0.29 — 
23.6) (0.40) vs Asian-Brazilian 
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Mixed 


Caucasian 


Asian-Brazilian 


Type 2 diabetes 


Yes 


No 


Hypertension 


Yes 


No 


Other 
respiratory 
diseases 


Yes 


No 


Cardiovascular 
diseases 


Yes 


No 


Cancer 


Yes 


No 


Smoking 


Yes 


30/1,113 
(2.7%) 


102/3,515 
(2.9%) 


1/90 
(1.1%) 


17/112 
(15.2%) 


120/4,747 
(2.5%) 


36/285 
(12.6%) 


101/4,574 
(2.2%) 


1/12 
(8.3%) 


136/4,837 
(2.8%) 


1/8 
(12.5%) 


136/4,851 
(2.8%) 


1/9 
(11.1%) 


136/4,850 
(2.8%) 


5/33 
(15.2%) 


132/4,836 
(2.7%) 


1/20 
(5.0%) 


136/4,839 
(2.8%) 


1/73 
(1.4%) 


136/4,786 


19/682 (2.8%) 


76/2,192 (3.5%) 


0/60 (0.0%) 


14/63 (22.2%) 


85/2,971 (2.9%) 


27/166 (16.3%) 


72/2,868 (2.5%) 


0/6 (0.0%) 


99/3,028 (3.3%) 


1/6 (15.7%) 


98/3,028 (3.2%) 


1/5 (20.0%) 


98/3,029 (3.2%) 


4/15 (26.7%) 


95/3.019 
(31.4%) 


1/12 (8.3%) 


98/3,022 (3.2%) 


1/47 (2.1%) 


98/2,987 (3.3%) 
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0/58 (0.0%) 


0/221 (0.0%) 


0/0 (0.0%) 


0/9 (0.0%) 


0/274 (0.0%) 


0/23 (0.0%) 


0/260 (0.0%) 


0/0 (0.0%) 


0/283 (0.0%) 


0/1 (0.0%) 


0/282 (0.0%) 


0/1 (0.0%) 


0/282 (0.0%) 


0/2 (0.0%) 


0/281 (0.0%) 


0/2 (0.0%) 


0/281 (0.0%) 


0/3 (0.0%) 


0/280 (0.0%) 


11/373 (2.9%) 


26/1,102 (2.4%) 


1/30 (3.3%) 


3/40 (7.5%) 


35/1,502 (2.3%) 


9/96 (9.4%) 


29/1,446 (2.0%) 


1/6 (16.7%) 


37/1,536 (2.4%) 


0/1 (0.0%) 


38/1,541 (2.5%) 


0/3 (0.0%) 


38/1,539 (2.5%) 


1/16 (6.3%) 


37/1,526 (24.2%) 


0/6 (0.0%) 


38/1,536 (2.5%) 


0/23 (0.0%) 


38/1,519 (%) 


0.14 


<0.0001 


n/a 


0.006 


<0.0001 


0.0003 


<0.0001 


n/a 


<0.0001 


n/a 


<0.0001 


n/a 


<0.0001 


0.058 


<0.0001 


n/a 


<0.0001 


n/a 


<0.0001 


0.93 (0.61 — 1.40) (0.72) vs Caucasians / 2.47 
(0.33 — 18.3) (0.38) vs Asian-Brazilian 


2.66 (0.37 — 19.3) (0.33) vs Asian-Brazilian 


6.90 (3.99 — 11.9) (<0.0001) 


6.40 (4.29 — 9.56) (<0.0001) 


3.14 (0.40 — 24.5) (0.27) 


4.95 (0.61 — 40.5) (0.14) 


4.33 (0.54 — 34.9) (0.17) 


6.36 (2.42 — 16.7) (0.0002) 


1.82 (0.24 — 13.7) (0.56) 


0.47 (0.07 — 3.44) (0.46) 
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(2.8%) 
History of 
9.56 (2.64 — 34.8) (0.0006) 
stroke 
3/14 
Yes 3/10 (30.0%) 0/1 (0.0%) 0/3 (0.0%) n/a 
(21.4%) 
134/4,845 
No (2.8%) 96/3,024 (3.2%) 0/282 (0.0%) 38/1,539 (2.5%) <0.0001 
. ‘0 
History of MI 26.4 (5.85 — 119.2) (<0.0001) 
Sm 
Yes 2/4 (50.0%) 0/0 (0.0%) 1/3 (33.3%) 0.087 
(42.9%) 
134/4,852 
No (2.8%) 97/3,030 (3.2%) 0/283 (0.0%) 37/1,539 (2.4%) <0.0001 
. ‘0 


TABLE 9: Pre-matching hospitalization rates in overall and in subpopulations among ivermectin 
non-users, irregular users, and regular users 


MI = myocardial infarction; COPD = chronic obstructive pulmonary disease; y/o = years old; Cl = confidence interval; n/a = not applicable. 


Hospitalization 


mine Non-users versus regular users Non-users versus irregular users Regular versus irregular users 
Unadjusted Multivariate Unadjusted Multivariate Unadjusted Multivariate 
PONNE hospital risk adjusted hospital hospital risk adjusted hospital hospital risk adjusted hospital 
3 ratio (95% CI) risk ratio (95% CI) ratio (95% CI) risk ratio (95% CI) ratio (95% CI) risk ratio (95% CI) 
and p-value (p) and p-value (p) and p-value (p) and p-value (p) and p-value (p) and p-value (p) 
Overal 0.05 (0.003 — 0.02 (0.00-0.27) 0.75 (0.51 — 0.68 (0.48-—0.96) 0.07 (0.004 — 0.03 (0.00 — 0.90) 
0.84) (0.037) (< 0.0001) 1.09) (0.14) (0.03) 1.13) (0.06) (0.043) 
Age 
aoia 7.12 (0.29 — 0.53 (0.00 — 844.1) 2.13 (0.13 — 1.73 (0.14 — 21.2) 3.35 (0.13 — 0.35 (0.00 — 
y 177.5) (0.23) (0.87) 34.1) (0.59) (0.67) 83.5) (0.46) 3625.6) (0.82) 
30-50 vio 0.24 (0.01 — 0.17 (0.01 - 5.00) 0.59 (0.25 — 0.58 (0.25 — 1.33) 0.40 (0.02 — 0.14 (0.001 — 14.3) 
% 4.00) (0.32) (0.30) 1.38) (0.22) (0.20) 7.03) (0.53) (0.40) 
ss0ylð 0.03 (0.002 — 0.03 (0.001 — 0.62) 0.64 (0.41 — 0.68 (0.46 — 1.01) 0.05 (0.003 — 0.02 (0.00 — 2.11) 
Z 0.51) (0.015) (0.024) 0.99) (0.047) (0.059) 0.81) (0.035) (0.099) 
Sex 
Female 0.13 (0.008 — 0.02 (0.00 — 4.20) 0.77 (0.43 — 0.69 (0.41 — 1.16) 0.17 (0.01 — 0.05 (0.00 — 3.60) 
2.15) (0.15) (0.15) 1.36) (0.36) (0.16) 2.82) (0.22) (0.17) 
Rais 0.08 (0.01 — 0.02 (0.00 — 1.86) 0.75 (0.45 — 0.67 (0.42 — 1.07) 0.11 (0.007 — 0.04 (0.00 — 2.37) 
1.34) (0.08) (0.088) 1.24) (0.26) (0.095) 1.81) (0.12) (0.12) 
Race 
Ean 2.38 (0.11 - als 0.29 (0.01 — 0.04 (0.00 — 32.4) 8.33 (0.15 — HE 
51.4) (0.58) 5.44) (0.40) (0.35) 473.6) (0.30) 
; 0.29 (0.02 — 1.06 (0.50 — 0.81 (0.40-0.98) 0.27 (0.02 — 
Mixed n/a n/a 
4.88) (0.39) 2.25) (0.88) (0.038) 4.63) (0.37) 
- 0.06 (0.004 — 0.67 (0.43 — 0.64 (0.43 - 0.98) 0.09 (0.01 — 
Caucasian n/a n/a 
1.01) (0.051) 1.06) (0.085) (0.038) 1.51) (0.095) 
6.15 (0.24 — 
Asian-Brazilian n/a n/a ( n/a (0.42) n/a n/a 
155.7) (0.27) 
Type 2 
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Cureus 


diabetes 

Yes 0.18 (0.01 — 
3.28) (0.25) 

No 0.06 (0.004 — 
0.99) (0.049) 

Hypertension 

vee 0.11 (0.01 — 
1.83) (0.12) 

No 0.07 (0.005 — 
1.20) (0.067) 

Asthma 

Yes n/a 

No 0.05 (0.003 — 
0.83) (0.037) 

COPD 

vee 1.22 (0.03 - 
48.2) (0.91) 

No 0.05 (0.003 — 
0.85) (0.038) 

Other 

respiratory 

diseases 

veš 1.00 (0.02 — 
40.3) (1.00) 

No 0.05 (0.003 — 
0.85) (0.038) 

Cardiovascular 

diseases 

ves 0.51 (0.02 — 
12.9) (0.68) 
0.05 (0.003 — 

No 
0.87) (0.04) 

Cancer 

vas 1.53 (0.05 — 
49.8) (0.81) 
0.05 (0.003 — 

No 
0.85) (0.038) 

Smoking 

vee 4.43 (0.15 — 
130.1) (0.39) 
0.05 (0.003 — 

No 
0.84) (0.038) 

History of 

stroke 

see 0.71 (0.02 — 
22.3) (0.85) 
0.05 (0.003 — 

No 
0.87) (0.039) 


0.01 (n/a) (0.58) 


0.04 (0.002 — 0.54) 
(0.016) 


0.01 (n/a) (0.36) 


0.04 (0.002 — 0.66) 
(0.025) 


n/a 


n/a 


0.28 (n/a) (0.97) 


0.02 (0.00 — 0.69) 
(0.031) 


0.01 (n/a) (0.67) 


0.03 (0.001 — 0.54) 
(0.018) 


0.01 (n/a) (0.82) 


0.03 (0.002 — 0.50) 
(0.015) 


0.04 (n/a) (0.92) 


0.01 (0.00 — 0.92) 
(0.046) 


0.08 (n/a) (0.89) 


0.03 (0.002 — 0.49) 
(0.014) 


0.02 (n/a) (0.82) 


0.02 (0.00 — 0.58) 
(0.022) 
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0.28 (0.08 — 
1.06) (0.061) 


0.34 (0.11 - 
1.04) (0.059) 


0.53 (0.24 — 
1.19) (0.12) 


0.79 (0.51 — 
1.23) (0.30) 


3.55 (0.12 — 
105.8) (0.47) 


0.73 (0.50 — 
1.07) (0.11) 


1.22 (0.03 - 
48.2) (0.91) 


0.75 (0.51 — 
1.09) (0.13) 


0.50 (0.03 — 
9.46) (0.64) 


0.76 (0.52 — 
1.11) (0.15) 


0.18 (0.02 — 
1.88) (0.15) 


0.76 (0.52 — 
1.12) (0.17) 


0.59 (0.02 — 
16.7) (0.76) 


0.76 (0.52 — 
1.11) (0.15) 


0.66 (0.03 — 
16.8) (0.80) 


0.76 (0.52 — 
1.11) (0.15) 


0.31 (0.01 — 
7.69) (0.47) 


0.77 (0.53 — 
1.13) (0.18) 


0.46 (0.14 — 1.49) 
(0.19) 


0.71 (0.49 — 1.03) 
(0.067) 


0.57 (0.28 — 1.15) 
(0.12) 


0.72 (0.48 — 1.09) 
(0.12) 


n/a (0.55) 


0.66 (0.46 — 0.94) 
(0.021) 


6.69 (0.31 — 145.6) 
(0.23) 


0.64 (0.45 — 0.92) 
(0.016) 


0.27 (0.00 — 759.0) 
(0.75) 


0.67 (0.47 — 0.96) 
(0.029) 


0.11 (0.01 — 1.67) 
(0.11) 


0.70 (0.49 — 0.99) 
(0.048) 


0.01 (n/a) (0.72) 


0.68 (0.48 — 0.97) 
(0.032) 


0.092 (0.00 — 
338.6) (0.57) 


0.68 (0.48 — 0.97) 
(0.035) 


n/a (0.73) 


0.69 (0.48 — 0.98) 
(0.038) 


0.56 (0.03 — 
11.9) (0.71) 


0.08 (0.005 — 
1.24) (0.07) 


0.20 (0.01 — 
3.49) (0.27) 


0.09 (0.01 — 
1.51) (0.095) 


n/a 


0.07 (0.004 — 
1.15) (0.063) 


1.00 (0.01 = 
92.4) (1.00) 


0.07 (0.004 — 
1.13) (0.061) 


2.33 (0.03 — 
182.9) (0.70) 


0.07 (0.004 — 
1.13) (0.061) 


2.07 (0.06 — 
66.3) (0.68) 


0.07 (0.004 — 
1.15) (0.063) 


2.60 (0.04 — 
170.4) (0.65) 


0.07 (0.004 — 
1.13) (0.061) 


6.71 (0.11 - 
396.2) (0.36) 


0.07 (0.004 — 
1.12) (0.06) 


2.33 (0.03 — 
182.9) (0.70) 


0.07 (0.004 — 
1.13) (0.06) 


0.01 (n/a) (0.62) 


0.05 (0.002 — 0.96) 
(0.047) 


0.02 (n/a) (0.44) 


0.04 (0.002 — 1.15) 
(0.061) 


n/a 


n/a 


0.01 (n/a) (0.87) 


0.04 (0.002 — 0.86) 
(0.04) 


n/a (0.99) 


0.04 (0.002 — 0.88) 
(0.041) 


0.38 (n/a) (0.97) 


0.04 (0.002 — 0.84) 
(0.038) 


n/a (0.95) 


0.04 (0.002 — 0.84) 
(0.038) 


n/a (0.99) 


0.04 (0.002 — 0.83) 
(0.037) 


n/a (0.97) 


0.04 (0.002 — 0.84) 
(0.038) 
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Cureus 


History of MI 
Yes n/a 

0.05 (0.003 — 
No 

0.86) (0.039) 


TABLE 10: Head-to-head comparisons of hospitalization rates between ivermectin non-users and 


n/a 


n/a 


0.50 (0.02 — 
11.1) (0.66) 


0.74 (0.51 — 
1.09) (0.13) 


0.17 (0.01 — 5.14) 
(0.31) 


n/a n/a 


0.69 (0.48-0.97) 0.07 (0.004 — 


(0.035) 4.15) (0.063) 2 


regular users, non-users and irregular users, and irregular and regular users 


MI = myocardial infarction; COPD = chronic obstructive pulmonary disease; y/o = years old; Cl = confidence interval; n/a = not applicable. 
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Cureus 


Citywide program of COVID-19 prophylaxis with 
BEFORE MATCHING . P sos f 
ivermectin (Itajaí, Brazil) 
NON REGULAR X IRREGULAR 
COVID-19 IVERMECTIN USE IVERMECTIN USE IVERMECTIN USE 
hospitalization 
rate 3,034 283 1,542 
99 hospitalizations O hospitalizations 38 hospitalizations 
Ov Il 3.3% 0.0% 2.5% 
era hospitalization rate hospitalization rate hospitalization rate 
population 
198% vs non-users 32% vs non-users 
(RR, 0.02; 95%C1, 0.00 - 0.27; p < 0.0001) (RR, 0,68; 95%CI, 0.48 — 0.96; p = 0.03) 
497% vs irregular users 
(RR, 0.03; 95%C!, 0.00 — 0.90; p = 0.043) 
42 hospitalizations O hospitalizations 17 hospitalizations 
(out of 1,624) (out of 141) (out of 853) 
2.6% 0.0% 2.0% 
hospitalization rate hospitalization rate hospitalization rate 
Females 
1.98% vs non-users 131% vs non-users 
(RR, 0.13; 95%CI, 0.00 — 4.20; p = 0.15) (RR, 0.69; 95%C1, 0,41 — 1.16; p = 0.16) 
495% vs irregular users 
(RR, 005; 95%, 000- 3.60; p= 0.17) 
57 hospitalizations O hospitalizations 21 hospitalizations 
{out of 1,410) (out of 142) (out of 689) 
4.0% 0.0% 3.0% 
hospitalization rate hospitalization rate hospitalization rate 
Males 
98% vs non-users 133% vs non-users 
IRR, 0.02; 95%, 0.00- 1.86; p = 0.088) IRR, 0,67; 951, 0.42 — 1.07; p = 0.095) 
1.96% vs irregular users 
IRR, 0.04; 95X1, 0.00- 2.37; p= 0.12) 
74 hospitalizations O hospitalizations 30 hospitalizations 
(out of 608) (out of 113) (out of 370) 
12.7% 0.0% 8.1% 
hospitalization rate hospitalization rate hospitalization rate 
Above 50 y/o 
197% vs non-users 32% vs non-users 
(RR, 0.03; 95%C1, 0.00 - 0.62; p = 0.024) (RR, 0.68; 95%CI, 0.46 — 1.01; p = 0.059) 
98% vs irregular users 
(RR, 0.02; 95%C1, 0.00- 2.11; p = 0,099) 
14 hospitalizations 0 hospitalizations 3 hospitalizations 
{out of 63) (out of 9) (out of 40) 
22.2% 0.0% 7.5% 
hospitalization rate | hospitalization rate hospitalization rate 
Type 2 
diabetes 99% vs non-users 54% vs non-users 
(RR, 0.01; 95%CI, n/a; p = 0.58) (RR, 0.46; 95%C1, 0.14 — 1.49; p = 0.19) 
499% vs irregular users 
(RR, 0.01; 95%CI, n/a; p = 0.28) 
27 hospitalizations O hospitalizations 9 hospitalizations 
(out of 166) {out of 23) (out of 96) 
16.3% 0.0% 9.4% 
hospitalization rate hospitalization rate hospitalization rate 
Hypertension 
YP 99% vs non-users 53% vs non-users 
(RR, 0.01; 95%C;, n/a; p=0.36) (RR, 0.47; 95%), 0.28~1.15; p= 0.12) 
498% vs irregular users 
IRR, 0.02; 951, n/a: p = 0.84) 
04 hospitalizations 0 hospitalizations 01 hospitalizations 
Cardiovascular (out of 15) (outot 2) | (outof16) | 
diseases 26.7% 0.0% 6.3% 
hospitalization rate hospitalization rate hospitalization rate 


Results are provided after multivariate adjusted analysis. 


FIGURE 7: Hospitalization rates before matching in ivermectin non- 
users, irregular users, and regular users 


RR = risk ratio; CI = confidence interval. 


Appendix 1B: mortality rates in non-users, irregular users, and regular 
users before matching 


Tables 11, 12 and Figure 8 describe mortality rates before matching in ivermectin non-users, irregular users, 
and regular users. Mortality rate was 2.6% in ivermectin non-users (79/3,034 deaths), 1.9% in irregular users 
(29/1,542 deaths), and 0.7% in regular users (2/283 deaths). Regular use of ivermectin was associated with 
85% reduction in mortality rate when compared to non-use of ivermectin (RR: 0.15; 95% CI: 0.04 - 0.59; p = 
0.007) and 77% reduction when compared to irregular users (RR: 0.23; 95% CI: 0.06 - 0.94; p = 0.031). 
Irregular use of ivermectin was associated with 36% reduction in mortality rate compared to non-use (RR: 
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Cureus 


Mortality 
rates 


Population 


Overall 
Age 
<30 y/o 


30-50 y/o 


>50 y/o 
Sex 


Female 


Male 


Race 


Afro-Brazilian 


Mixed 


Caucasian 


Asian-Brazilian 


Type 2 
diabetes 


Yes 


No 
Hypertension 


Yes 


No 
Asthma 


Yes 


COPD 


Non-matched (pre-matching) groups 


Overall (n 
= 4,859) 


110/4,859 
(2.3%) 


1/1,280 
(0.1%) 


12/2,488 
(1.2%) 


97/1,091 
(8.9%) 


51/2,618 
(1.9%) 


59/2,241 
(2.6%) 


41141 
(2.8%) 


19/1,113 
(1.7%) 


85/3,515 
(2.4%) 


2/90 
(2.2%) 


20/112 
(17.9%) 


90/4,747 
(1.9%) 


36/285 
(12.6%) 


74/4,574 
(1.6%) 


2/12 
(16.7%) 


108/4,837 
(2.2%) 


0.64; 95% CI: 0.43 - 0.96; p = 0.041). 


Non-users - 
deaths (%) 
(n = 3,034) 


79/3,034 
(2.6%) 


1/844 
(0.1%) 


10/1,582 
(0.6%) 


68/608 
(11.1%) 


36/1,624 
(2.2%) 


43/1,410 
(3.1%) 


4/100 
(4.0%) 


12/682 
(1.8%) 
62/2,192 


(2.8%) 


1/60 (1.7%) 


16/63 
(25.4%) 


63/2,971 
(2.1%) 


28/166 
(16.9%) 


51/2,868 
(1.8%) 


1/6 (16.7%) 


78/3,028 
(2.6%) 


Regular 
users - 
deaths (%) 
(n = 283) 


2/283 
(0.7%) 


0/39 (0.0%) 


0/131 
(0.0%) 


2/113 
(1.8%) 


0/141 
(0.0%) 


2/142 
(1.4%) 


0/4 (0.0%) 


0/58 (0.0%) 


2/221 
(0.9%) 


0/0 (0.0%) 


1/9 (11.1%) 


1/274 
(0.4%) 


1/23 (4.3%) 


1/260 
(0.4%) 


0/0 (0.0%) 


2/283 
(0.7%) 
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Irregular 
users - 
deaths (%) (n 
= 1,542) 


29/1,542 
(1.9%) 


0/397 (0.0%) 


2/775 (2.6%) 
27/370 


(7.3%) 


15/853 
(1.8%) 


14/689 
(2.0%) 


0/37 (0.0%) 


7/373 (1.9%) 


21/1,102 
(1.9%) 


1/30 (3.3%) 


3/40 (7.5%) 


26/1,502 
(1.7%) 


7/96 (7.3%) 


22/1 ,446 
(1.5%) 


1/6 (16.7%) 


28/1,536 
(18.2%) 


Overall p-values 


Unadjusted 
overall p-value 
(between three 


groups) 


0.058 


0.05 


0.002 


<0.0001 


0.25 


0.001 


0.4 


0.085 


0.31 


0.059 


0.11 


0.037 


0.22 


0.91 


0.054 


Unadjusted overall p-value (p) and risk ratio 
(95% Cl) between subpopulations 


<0.0001 


6.20 (0.81 — 47.7) (0.08) vs < 30 y/o 


13.8 (1.92 — 99.0) (0.009) vs < 30 y/o / 2.25 
(1.22 — 4.06) (0.009) vs 30-50 y/o 


1.36 (0.93 — 1.99) (0.11) vs female 


0.54 


1.18 (0.43 — 3.26) (0.75) vs Caucasians / 
1.68 (0.56 — 5.01) (0.35) vs mixed / 1.28 
(0.23 — 7.16) (0.78) vs Asian-Brazilian 


0.70 (0.42 — 1.16) (0.17) vs Caucasians / 
0.76 (0.18 — 3.33) (0.72) vs Asian-Brazilian 


1.09 (0.26 — 4.50) (0.90) vs Asian-Brazilian 


11.2 (6.64 — 19.0) (<0.0001) 


8.79 (5.79 — 13.4) (<0.0001) 


8.76 (1.90 — 40.4) (0.005) 


6.21 (0.76 — 51.0) (0.088) 
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Cureus 


Yes 


No 


Other 
respiratory 
diseases 


Yes 


No 


Cardiovascular 
diseases 


Yes 


No 
Smoking 


Yes 


No 


History of 
stroke 


Yes 


No 
History of MI 


Yes 


1/8 
1/6 (16.7%) 0/1 (0.09 
(12.5%) is er) 
109/4,851 78/3,028 2/282 
(2.2%) (2.6%) (0.7%) 
13 1/5 (20.0%) 0/1 (0.0%) 
(11.1%) A 7 
109/4,850 77/3,029 2/282 
(2.2%) (2.6%) (0.7%) 
5/33 2/115 
0/2 (0.0%) 


(15.2%) (13.3%) 


105/4,836 77/3,019 2/281 


(2.2%) (2.6%) (0.7%) 
1/20 1/12 (8.3%) 0/2 (%) 
(5.0%) Ti 7 
109/4,839 78/3,022 2/281 
(2.2%) (2.6%) (0.7%) 
es 1/47 (2.1%) 0/3 (0.0%) 
(1.4%) -170 -0% 
109/4,786 78/2,987 2/280 
(2.2%) (2.6%) (0.7%) 
3/14 4/10 


0/1 (0.0% 
(21.4%) (40.0%) eo”) 
107/4,845 75/3024 2/282 
(2.2%) (2.5%) (0.7%) 


3/7 


1/4 (25.0%) 0/0 (0.09 
(42.9%) 0%) (0.020) 


107/4,852 78/3,030 2/283 
(2.2%) (2.6%) (0.7%) 


0/1 (0.0%) 


29/1,541 
(1.9%) 


0/3 (0.0%) 


29/1,539 
(1.9%) 


0/16 (0.0%) 


29/1,526 
(1.9%) 


0/6 (0.0%) 


29/1,536 
(1.9%) 


0/23 (0.0%) 


29/1,519 
(1.9%) 


0/3 (0.0%) 


29/1,539 
(1.9%) 


0/3 (0.0%) 


29/1,539 
(1.9%) 


n/a 


0.065 


5.44 (0.67 — 43.8) (0.11) 


n/a 


0.074 


8.05 (3.05 — 21.2) (<0.0001) 


0.001 


0.078 


2.28 (0.30 — 17.2) (0.42) 


0.032 


0.066 


0.59 (0.08 — 4.33) (0.61) 


0.89 


0.064 


12.1 (3.32 — 43.9) (0.0002) 


0.001 


0.093 


33.3 (7.35 — 150.4) (<0.0001) 


n/a 


0.065 


TABLE 11: Pre-matching mortality rates in overall and in subpopulations among ivermectin non- 


users, irregular users, and regular users 


MI = myocardial infarction; COPD = chronic obstructive pulmonary disease; y/o = years old; CI = confidence interval; n/a = not applicable. 


Mortality 
rates 


Non-users versus regular users 


Unadjusted Multivariate 
mortality risk adjusted mortality 
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Regular users versus irregular 


Non-users versus irregular users 


Unadjusted 
mortality risk 


users 
Multivariate Unadjusted Multivariate 
adjusted mortality mortality risk adjusted mortality 
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Cureus 


Population 


Overall 
Age 


<30 y/o 


30-50 y/o 


>50 y/o 
Sex 


Female 


Male 
Race 


Afro-Brazilian 


Mixed 


Caucasian 


Asian-Brazilian 


Type 2 
diabetes 


Yes 


No 
Hypertension 


Yes 


No 
Asthma 


Yes 


No 


Other 
respiratory 
diseases 


ratio (95% Cl) 
and p-value (p) 


0.27 (0.07 - 
1.09) (0.066) 


7.11 (0.29 - 
177.5) (0.23) 


0.57 (0.03 — 
9.77) (0.70) 


0.14 (0.03 — 
0.59) (0.007) 


0.15 (0.01 — 
2.52) (0.19) 


0.45 (0.11 — 
1.89) (0.28) 


2.38 (0.11 - 
51.4) (0.58) 


0.46 (0.03 — 
7.84) (0.59) 


0.31 (0.08 — 
1.29) (0.11) 


n/a 


0.37 (0.04 — 
3.17) (0.36) 


0.17 (0.02 - 
1.22) (0.078) 


0.22 (0.03 — 
1.73) (0.15) 


0.21 (0.03 — 
1.55) (0.13) 


3.67 (0.05 — 
274.5) (0.56) 


0.27 (0.07 - 
1.10) (0.068) 


1.22 (0.03 - 
48.2) (0.91) 


0.27 (0.07 - 
1.11) (0.069) 


risk ratio (95% Cl) 
and p-value (p) 


0.15 (0.04 — 0.59) 
(0.007) 


1.02 (0.00 - 
1069.7) (0.99) 


0.33 (0.01 — 13.5) 
(0.56) 


0.15 (0.04 — 0.64) 
(0.01) 


0.08 (0.004 — 1.47) 
(0.089) 


0.20 (0.04 — 0.99) 
(0.048) 


n/a 


n/a 


n/a 


n/a 


0.46 (0.07 — 3.27) 
(0.44) 


0.12 (0.02 — 0.63) 
(0.012) 


0.24 (0.03 — 1.71) 
(0.15) 


0.11 (0.01 — 0.76) 
(0.025) 


n/a 


/a 


=l 


0.92 (n/a) (1.00) 


0.15 (0.04 — 0.61) 
(0.007) 
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ratio (95% CI) 
and p-value (p) 


0.72 (0.47 — 
1.10) (0.13) 


0.71 (0.03 — 
17.4) (0.83) 


0.41 (0.09 — 
1.86) (0.25) 


0.63 (0.39 — 
0.997) (0.048) 


0.79 (0.43 — 
1.45) (0.45) 


0.66 (0.36 — 
1.21) (0.18) 


0.29 (0.02 — 
5.44) (0.41) 


1.07 (0.42 = 
2.74) (0.89) 


0.67 (0.40 — 
1.11) (0.11) 


2.03 (0.13 — 
33.7) (0.62) 


0.24 (0.06 — 
0.88) (0.031) 


0.81 (0.51 — 
1.29) (0.38) 


0.39 (0.16 — 
0.93) (0.033) 


0.85 (0.52 — 
1.41) (0.54) 


1.00 (0.05 - 
20.8) (1.00) 


0.70 (0.45 — 
1.09) (0.11) 


1.22 (0.03 - 
48.2) (0.91) 


0.73 (0.47 — 
1.12) (0.14) 


risk ratio (95% Cl) 
and p-value (p) 


0.64 (0.43 — 0.96) 
(0.041) 


0.49 (0.01 — 18.4) 
(0.70) 


0.47 (0.12 — 1.84) 
(0.28) 


0.68 (0.45 — 1.03) 
(0.071) 


0.70 (0.40 — 1.24) 
(0.22) 


0.59 (0.34 — 1.03) 
(0.065) 


0.05 (0.00 — 36.0) 
(0.36) 


0.76 (0.31 — 1.86) 
(0.55) 


0.64 (0.40 — 1.01) 
(0.058) 


1.81 (0.13 — 24.5) 
(0.65) 


0.42 (0.14 — 1.33) 
(0.14) 


0.70 (0.45 — 1.08) 
(0.11) 


0.48 (0.22 — 1.01) 
(0.054) 


0.75 (0.46 — 1.20) 
(0.23) 


4.34 (0.04 — 446.1) 
(0.53) 


0.63 (0.42 — 0.94) 
(0.024) 


4.98 (0.34 — 73.9) 
(0.24) 


0.63 (0.34 — 094) 
(0.024) 


ratio (95% Cl) 
and p-value (p) 


0.36 (0.09 — 
1.51) (0.16) 


10.2 (0.20 - 
505.9) (0.24) 


1.21 (0.06 - 
25.1) (0.90) 


0.24 (0.06 — 
0.98) (0.046) 


0.19 (0.01 — 
3.13) (0.24) 


0.67 (0.15 — 
2.89) (0.59) 


4.80 (0.11 - 
214.5) (0.42) 


0.42 (0.02 — 
7.31) (0.55) 


0.45 (0.11 — 
1.89) (0.28) 


n/a 


1.56 (0.18 — 
13.3) (0.69) 


0.20 (0.03 — 
1.48) (0.12) 


0.54 (0.07 — 
4.13) (0.56) 


0.25 (0.03 — 
1.84) (0.17) 


2.67 (0.23 — 
30.4) (0.43) 


0.37 (0.09 — 
1.56) (0.18) 


0.50 (0.04 — 
7.10) (0.61) 


0.38 (0.09 — 
1.57) (0.18) 


risk ratio (95% Cl) 
and p-value (p) 


0.23 (0.06 — 0.94) 
(0.031) 


1.30 (0.00 = 
4060.2) (0.95) 


0.53 (0.01 — 34.1) 
(0.76) 


0.24 (0.04 — 1.05) 
(0.058) 


0.10 (0.004 — 2.59) 
(0.17) 


0.42 (0.10 — 1.86) 
(0.26) 


n/a 


n/a 


n/a 


n/a 


0.90 (0.11 — 7.59) 
(0.92) 


0.17 (0.03 — 1.03) 
(0.054) 


0.53 (0.07 — 4.29) 
(0.95) 


0.20 (0.03 — 1.19) 
(0.077) 


n/a 


n/a 


0.03 (n/a) (0.89) 


0.27 (0.07 — 1.04) 
(0.058) 
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= 1.00 (0.02 — 
40.3) (1.00) 
0.27 (0.07 — 
No 


1.11) (0.069) 


Cardiovascular 


diseases 
ven 1.08 (0.04 — 
30.0) (0.96) 
0.27 (0.07 — 
No 
1.12) (0.072) 
Cancer 
vaš 1.53 (0.05 — 
49.8) (0.81) 
0.27 (0.07 — 
No 
1.11) (0.069) 
Smoking 
ine 4.43 (0.15 — 
130.1) (0.39) 
0.27 (0.07 — 
No 
1.10) (0.067) 
History of 
stroke 
ae 0.48 (0.02 — 
14.7) (0.68) 
0.28 (0.07 — 
No 
1.15) (0.077) 
History of MI 
Vee 0.78 (0.02 — 
32.4) (0.89) 
0.27 (0.07 — 
No ( 


1.10) (0.068) 


0.04 (n/a) (0.48) 


0.17 (0.05 — 0.62) 
(0.01) 


0.03 (n/a) (0.80) 


0.16 (0.04 — 0.62) 
(0.008) 


0.3 (n/a) (0.95) 


0.16 (0.04 — 0.61) 
(0.008) 


0.33 (n/a) (0.92) 


0.15 (0.04 — 0.61) 
(0.008) 


0.005 (0.00 - > 
1000) (0.64) 


0.20 (0.06 — 0.70) 
(0.012) 


n/a 


n/a 


0.43 (0.01 — 
14.1) (0.63) 


0.74 (0.48 — 
1.13) (0.16) 


0.16 (0.01 — 
3.71) (0.26) 


0.74 (0.48 — 
1.14) (0.17) 


0.59 (0.02 — 
16.7) (0.76) 


0.73 (0.47 — 
1.12) (0.15) 


0.66 (0.03 — 
16.8) (0.80) 


0.73 (0.47 — 
1.12) (0.15) 


0.21 (0.01 — 
5.05) (0.33) 


0.76 (0.49 — 
1.16) (0.20) 


0.33 (0.01 — 
11.3) (0.54) 


0.73 (0.47 — 
1.12) (0.15) 


0.13 (n/a) (0.78) 


0.65 (0.44 — 0.87) 
(0.037) 


0.02 (n/a) (0.41) 


0.69 (0.46 — 1.02) 
(0.061) 


0.01 (n/a) (0.74) 
0.66 (0.44 — 0.98) 


(0.039) 


0.12 (0.00 — 272.7) 
(0.59) 


0.66 (0.44 — 0.98) 
(0.04) 


0.002 (n/a) (0.71) 


0.68 (0.46 — 1.02) 
(0.065) 


0.002 (n/a) (0.74) 


0.67 (0.45 — 0.99) 
(0.047) 


2.0 (0.06 — 68.5) 
(0.70) 


0.36 (0.09 — 
1.51) (0.16) 


5.67 (0.14 — 
234.5) (0.36) 


0.36 (0.09 — 
1.50) (0.16) 


2.33 (0.06 — 
92.4) (0.65) 


0.36 (0.09 — 
1.51) (0.16) 


6.25 (0.14 — 
272.5) (0.23) 


0.36 (0.09 — 
1.50) (0.16) 


2.00 (0.06 — 
68.5) (0.70) 


0.36 (0.09 — 
1.51) (0.16) 


4.00 (0.15 — 
103.2) (0.40) 


0.36 (0.09 — 
1.51) (0.16) 


0.26 (0.07 — 1.03) 
(0.054) 


0.98 (n/a) (1.00) 


0.61 (n/a) (0.98) 


0.26 (0.07 — 1.03) 
(0.056) 


14.7 (n/a) (0.92) 


0.26 (0.07 — 1.03) 
(0.055) 


1.69 (n/a) (0.97) 


0.27 (0.07 — 1.03) 
(0.056) 


1.73 (n/a) (0.99) 


0.26 (0.07 — 1.03) 
(0.055) 


TABLE 12: Head-to-head comparisons of mortality rates between ivermectin non-users and 
regular users, non-users and sporadic users, and between sporadic and regular users 


MI = myocardial infarction; COPD = chronic obstructive pulmonary disease; y/o = years old; Cl = confidence interval; n/a = not applicable. 
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rari Citywide program of COVID-19 prophylaxis with 
ivermectin (Itajaí, Brazil) 
X REGULAR X IRREGULAR 
COVID-19 NONIIVERMECTIN USE IVERMECTIN USE IVERMECTIN USE 
mortality 
rate 3,034 283 1,542 
J l | 
+ 
79 deaths 02 deaths | 29 deaths 
2.6% 0.7% | 1.9% 
Overall mortality rate mortality rate mortality rate 
population 85% vs non-users 136% vs non-users 
(RR, 0.15; 95%CI, 0.04 ~ 0,59; p = 0.007) (RR, 0.64; 95%CI, 0.43 - 0.96; p = 0.041) 
77% vs irregular users 
IRR, 023; 95%, 006 0.84; p = 0.031) 
36 deaths Odeaths 15 deaths 
(out of 1,624) (out of 141) (out of 853) 
2.2% 0.0% 1.8% 
mortality rate mortality rate mortality rate 
somons 492% vs non-users 130% vs non-users 
(RR, 0,08; 95%C1, 0,00 — 1.47; p = 0,089) (RR, 0.70; 95%CI, 0.40 — 1.24; p = 0.22) 
490% vs irregular users 
(RR, 0.10; 95%C1, 0.004 — 2.59; p = 0.17) 
43 deaths 02 deaths 14 deaths 
(out of 1,410) (out of 142) (out of 689) 
3.1% 1.4% 2.0% 
mortality rate mortality rate mortality rate 
Males 80% vs non-users 41% vs non-users 
(RR, 0.20; 95%CI, 0.04 — 0.99; p = 0.048) {RR, 0.59; 95%Cl, 0.34 — 1.03; p = 0.065) 
458% vs irregular users 
(RR, 0.42; 95%C!, 0.10 - 1,86; p = 0.26) 
68 deaths 02 deaths 27 deaths 
(out of 608) (out of 113) (out of 370) 
11.1% 1.8% 7.3% 
mortality rate mortality rate | mortality rate | 
Above 50 y/o 85% vs non-users 32% vs non-users 
(RR, 0.15; 95%CI, 0.04 ~ 0.64; p = 0.01) (RR, 0.68; 95%CI, 0.45 ~ 1.03; p = 0,071) 
176% vs irregular users 
(RR, 024; 95%CI, 0.04- 1.05; p = 0.058) 
16 deaths 01 death 03 deaths 
(out of 63) (out of 09) (out of 40) 
25.4% 11.1% 7.5% 
mortality rate mortality rate mortality rate 
Type 2 
diabetes 54% vs non-users 58% vs non-users 
(RR, 0.46; 95%CI, 0,10 — 4.19; p = 0.64) (RR, 0.42; 95%CI, 0.14 ~ 1.33; p = 0.14) 
410% vs irregular users 
(RR, 0.90;95%C1,0.11-7.59; p = 0.92) 
28 deaths 01 death 07 deaths 
(out of 166) (out of 23) (out of 96) 
16.9% 4.3% 7.3% 
mortality rate | mortality rate _ mortality rate 
Hypertension 
ypa 76% vs non-users 52% vs non-users 
(RR, 0.24; 95%CI, 0,03 ~ 1.71; p= 0.15) (RR, 0.48; 95%CI, 0,22 — 1.01; p = 0.054) 
47% vs irregular users 
ARR, 0.53; 95%CI, 0.07 - 4.29; p = 0.95) 
02 deaths Odeaths O deaths 
Cardiovaecul (out of 15) (out of 2) (out of 16) 
ar diseases Tak es oo% 
mortality rate mortality rate mortality rate 


Results are provided after multivariate adjusted analysis 


FIGURE 8: Pre-matching mortality rates in ivermectin non-users, 
irregular users, and regular users 


RR = risk ratio; CI = confidence interval. 


Additional Information 
Disclosures 


Human subjects: Consent was obtained or waived by all participants in this study. National Research Ethics 
Council (CONEP) issued approval 4.821.082. The study design, institutional review board (IRB) approval, and 
data analysis were done upon completion of the program. The study of the COVID-19 cases reported in the 
city of Itajaí (n = 9,956, including cases that occurred before July 7, 2020, as a comparison) was approved by 
the National Research Ethics Council (CONEP) (approval number: 4.821.082; protocol (CAAE) number: 
47124221.2.0000.5485). Animal subjects: All authors have confirmed that this study did not involve animal 
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subjects or tissue. Conflicts of interest: In compliance with the ICMJE uniform disclosure form, all authors 
declare the following: Payment/services info: Lucy Kerr is a paid consultant for Vitamedic, an ivermectin 
manufacturer, and Médicos Pela Vida (MPV), an organization that promotes ivermectin as a treatment for 
COVID-19. Flavio A. Cadegiani has performed paid consulting work for Vitamedic, an ivermectin 
manufacturer. Pierre Kory receives fees from COVID-19 Critical Care Alliance (FLCCC), an organization that 
promotes ivermectin as a treatment for COVID-19. Juan J. Chamie-Quintero is a paid contributor for the 
COVID-19 Critical Care Alliance (FLCCC), an organization that promotes ivermectin as a treatment for 
COVID-19. Financial relationships: Pierre Kory, Juan J. Chamie-Quintero declare(s) personal fees from 
Front Line COVID-19 Critical Care Alliance (FLCCC). Flavio A. Cadegiani, Pierre Kory declare(s) non- 
financial support from Front Line COVID-19 Critical Care Alliance (FLCCC). Lucy Kerr declare(s) non- 
financial support from Médicos Pela Vida (MPV). Jennifer A. Hibberd declare(s) non-financial support from 
Canadian Covid Care Alliance. Jennifer A. Hibberd declare(s) non-financial support from World Council for 
Health (WCH). Flavio A. Cadegiani, Lucy Kerr declare(s) personal fees from Vitamedic. Other relationships: 
Lucy Kerris is a paid consultant for both Vitamedic, an ivermectin manufacturer, and is co-founder, as well 
as acting as a paid consultant, for Médicos Pela Vida (MPV), an organization that promotes ivermectin as a 
treatment for COVID-19 and discourages COVID-19 vaccination. Flavio A. Cadegiani was a paid consultant 
(USD 1,600.00) for Vitamedic, an ivermectin manufacturer. Dr. Cadegiani is a founding member of the Front 
Line COVID-19 Critical Care Alliance (FLCCC), an organization that promotes ivermectin as a treatment for 
COVID-19. Pierre Kory is the President and Chief Medical Officer of the Front Line COVID-19 Critical Care 
Alliance (FLCCC), an organization that promotes ivermectin as a treatment for COVID-19 and discourages 
COVID-19 vaccination. Dr. Kory reports receiving payments from FLCCC. In February of 2022, Dr. Kory 
opened a private telehealth fee-based service to evaluate and treat patients with acute COVID, long-haul 
COVID, and post-vaccination syndromes. Dr. Kory expresses personal points of view against COVID-19 
vaccination. Jennifer A. Hibberd is a co-founder of the Canadian Covid Care Alliance and World Council for 
Health, both of which discourage vaccination and encourage ivermectin as a treatment for COVID-19. Juan J. 
Chamie-Quintero is a contributor to the Front Line COVID-19 Critical Care Alliance (FLCCC), lists the 
FLCCC as his employer on his LinkedIn page, and expresses personal points of view against COVID-19 
vaccination. 
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